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Abstract. A detailed analysis of transverse momentum spectra of several identified hadrons in high energy 
collisions within the canonical framework of the statistical model of hadronisation is performed. The study 
of particle momentum spectra requires an extension of the statistical model formalism used to handle 
particle multiplicities, which is described in detail starting from a microcanonical treatment of single 
hadronising clusters. Also, a new treatment of extra strangeness suppression is presented which is based 
on the enforcement of fixed numbers of ss pairs in the primary hadrons. The considered centre-of-mass 
energies range from ~ 10 to 30 GeV in hadronic collisions (np, pp and Kp) and from ~ 15 to 35 GeV in 
e + e~ collisions. The effect of the decay chain following hadron generation is accurately and exhaustively 
taken into account by a newly proposed numerical method. The exact p T conservation at low energy and 
the increasing hard parton emission at high energy bound the validity of the presently taken approach 
within a limited centre-of-mass energy range. However, within this region, a clear consistency is found 
between the temperature parameter extracted from the present analysis and that obtained from fits to 
average hadron multiplicities in the same collision systems. This finding indicates that in the hadronisation 
process the production of different particle species and their momentum spectra are two closely related 
phenomena governed by one parameter. 
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1 Introduction 

The idea of a statistical approach to hadron production 
in high energy collisions dates back to '50s Q and '60s 
|Q] and it has been recently revived by the observation 
that hadron multiplicities in e + e _ and pp collisions agree 
very well with a thermodynamical-like ansatz This 
finding has also been confirmed in hadronic collisions and 
it has been interpreted in terms of a pure statistical fill- 
ing of multi-hadronic phase space of assumed pre-hadronic 
clusters (or fireballs) formed in high-energy collisions, at a 
critical value of energy density J^,^,[^]. In this framework, 
temperature and other thermodynamical quantities have 
a purely statistical meaning and do not involve the exis- 
tence of a hadronic thermalisation process through multi- 
ple collisions on an event-by-event basis. Stated otherwise, 
statistical equilibrium is an intrinsic feature of the hadro- 
nisation process and hadrons are directly created in such 
a state ||]fj,|7|], as it was envisaged by Hagedorn ||. 

So far, this proposed statistical cluster hadronisation 
model has been mainly tested against measured abun- 
dances of different hadron species for a twofold reason. 
Firstly, unlike momentum spectra, they are quantities 
which are not affected by hard (perturbative) QCD dy- 



namical effects but are only determined by the hadronisa- 
tion process; indeed, in the framework of a multi-cluster 
model, they are Lorentz invariant quantities which are in- 
dependent of cluster's overall momentum. Secondly, they 
are fairly easy to calculate and provide a very sensitive 
test of the model yielding an accurate determination of 
the temperature. However, in order to establish the valid- 
ity of the model, it is necessary to test further observables 
and to assess their consistency with the results obtained 
for multiplicities. One of the best suited observables in this 
regard is the transverse momentum of identified hadrons 
(where transverse is meant to be with respect to beam line 
in high energy hadronic collision, and thrust or event axis 
in high energy e + e~ collisions) because, amongst all pro- 
jections of particle momentum, this is supposed to be the 
most sensitive to hadronisation or, conversely, the least 
sensitive to perturbative QCD dynamics. 

Actually, it has been known for a long time that trans- 
verse momentum spectra follow a Boltzmann distribu- 
tion in hadronic collisions and this very observation was 
pointed out by Hagedorn as a major indication in favour of 
statistical hadron production ||. It must be emphasized 
that the prediction of a thermal-like shape in principle 
only applies to particles directly emitted from the hadro- 
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nising source, whereas measured spectra also include par- 
ticles produced by decays of heavier hadrons. Yet, due 
to involved calculations into play, most analyses do not 
take into account the distortion of primordial hadronisa- 
tion spectrum induced by hadronic decays and try to fit 
the data straight through it. This problem has been dis- 
cussed in literature [ flO| and an analytical calculation has 
been developed to take into account the effect of two and 
three body decays |ll|,^2[ , which has then been used both 
for pp |12j and heavy ion collision Q[l|,[n|[l4j including 
most abundant resonances. In this paper we introduce a 
new method allowing to rigorously and exhaustively de- 
termine the contribution of all particle decays. Thence, by 
taking advantage of this technique, we have performed an 
analysis of many measured transverse momentum spec- 
tra of identified hadrons in a wide range of centre-of-mass 
energies for several kind of collisions. 

The paper is organised as follows: in Sect. 2 the sta- 
tistical hadronisation model is described in detail start- 
ing from a basic microcanonical treatment of clusters. In 
Sect. 3 the analytical formulae for transverse momentum 
spectra of hadrons within the statistical hadronisation 
model are derived whilst in Sect. 4 they are worked out for 
the primary and secondary component separately and the 
numerical method used to calculate them is described. In 
Sect. 5 and 6 the data analysis is presented and discussed. 
Conclusions are drawn in Sect. 7. 



2 Statistical hadronisation model 

The statistical hadronisation model assumes that in high 
energy collisions, as a consequence of strong interaction 
dynamics, a set of colour singlet clusters (or fireballs) en- 
dowed with mass, volume, internal quantum numbers and 
momentum, whose distribution is governed by the dy- 
namical stage of the process, is formed. These clusters 
are supposed to give rise to hadrons according to a pure 
statistical law in the multi-hadronic phase space defined 
by their four-momentum, volume and quantum numbers. 
This approach differs from another popular cluster hadro- 
nisation model mainly because clusters are provided 
with a volume, so that hadron production is governed 
by proper phase space rather than relativistic momentum 
space. In this framework, the use of statistical mechan- 
ics and thermodynamical quantities, such as temperature, 
which needs spacial dimension besides momentum space 
in order to be meaningful, is allowed. We emphasize once 
more that the introduction of such quantities does not en- 
tail any thermalisation process of hadrons after their for- 
mation and that statistical equilibrium simply means that 
all final available quantum mechanical states are equally 
likely. 

In refs. ||[l6| the statistical hadronisation model was 
described within a canonical framework, with clusters 
characterized by temperature and volume instead of mass 
and volume. Therein, it was shown that a particular 
choice of the probabilities of distributing quantum num- 
bers among the clusters and the assumption of a common 
temperature lead to a very simple expression for particle 



multiplicities and global particle correlations which are in 
fact the same as those relevant to one cluster having as 
volume the sum of the volumes of all clusters in their own 
rest frame. Furthermore, no dependence on single clus- 
ter properties nor on their number is left. The question 
arises whether such clusters are actually large enough to 
allow a canonical description. In fact, in principle, their 
mass and volume might be so small to require a more ap- 
propriate microcanonical framework, enforcing the exact 
conservation of energy and momentum in the calculation 
of the available multi-hadronic phase space. In the follow- 
ing, we will prove that a similar reduction of the expression 
of particle multiplicities holds in the microcanonical case, 
provided that suitable clusters configuration probabilities 
occur. 

2.1 From microcanonical to canonical ensemble 

In the canonical framework of the statistical hadronisa- 
tion model, the main tool to derive physical observables 
is the partition function Z, which is the sum over all 
physical states with fixed quantum numbers weighted by 
exp[— E/T], where E is the energy of the state and T the 
temperature []. Similarly in the microcanonical case, it is 
possible to introduce the sum over all physical states with 
fixed values of energy-momentum and quantum numbers, 
i.e. the density of states Q which, for the z th cluster, reads: 

= E - <*Q ,Q i = 1, ...,N (I) 

states 

where P, is the four-momentum of the cluster, Vi is its 
volume in the frame where four-momentum is Pi and 
Qi = (Qii) ■ ■ ■ , Qin) is a vector of its n relevant quantum 
numbers; Pi- S and Qi ;s are the corresponding quantities of 
a general multi-hadronic state which, in the ideal hadron 
gas approximation, is described by a set {rijk} of occu- 
pation numbers for each hadronic species j and for each 
phase space cell k, running from to I for fermions and 
from to oo for bosons. Hence: 

Pi;a = y^.Pjkrijk 

Jk 

Qi s = qjfijk (2) 

jk 

where is the quantum number vector for the j th hadron 
and pjk its four-momentum in the fc th phase space cell. 
The quantum numbers are supposed to be either integer- 
valued additive conserved quantities in strong interac- 
tions (namely electric charge, baryon number, strangeness, 

1 in a covariant formulation this weight is to be replaced by 
exp[— j3 ■ P] where P is the four-momentum and /3 the four- 
temperature vector. 

2 Throughout the paper, by {Ai} we mean a shorthand for 
the vector (Ai, . . . , An), either finite or infinite-dimensional. 
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charm and beauty) or positive integer-valued absolute 
numbers of valence quarks plus antiquarks. It is worth 
remarking that fli defined in Eq. ([!]) is a number of states 
per unit four-momentum cell, thus it is a Lorentz invari- 
ant quantity and has to depend on Lorentz scalars only. 
Otherwise stated: 



the function F finally reads: 
F{xi,<t>i) = 

-^3^(2J,- + 1) J d 3 p logtlie-'fe'^^^ 1 



y. 00 r 

QiiP^QiuVi) = C2 t (PiQ u V[) = Oi{P*MuV*) (3) = + *) E^ 1 )"^ d3 P e ~ 



where P/ is the Lorentz-transformed of Pi and V( is the 
volume in the corresponding frame, while P* and V* are 
the four-momentum (Mi,0) and the volume respectively 
in the cluster's rest frame. 

fli in Eq. ([!]) can be transformed by means of the in- 
tegral representations of Dirac's and Kronecker's delta: 



<5 4 (P - Pi-, s ) = J d 4 ^ exp [i(P< - P ;s ) ■ x t ] 

S QuQi,» = T^w J - J ex P MQi - Q ' 0;] 

(4) 

By plugging Eqs. (0) into Eq. (|j), the density of states 
reads: 



Qi = (2tt) 4 +" / d4xi J 6XP ' Xl + ^ ' ^ 



x II CX P [ - i n jk(p 3 k ■ Xi + q 3 ■ 4> i 

{n jk } jk 



(5) 



The calculation now proceeds by taking advantage of the 
commutability between sum and product in Eq. rtq) . How- 
ever, the sum over occupation number does not converge 
to a finite value for bosons as runs from to 00. The 
convergence is recovered if the time component of Xi is 
provided with a negative imaginary part —is. If we intro- 
duce such a term in Eq. (Q), then the bosonic sums can 
be performed and fli reads: 



fli = lim 



1 



dx° / d 3 x 



- (27r) 4+n 
x exp [i (Pi ■ xi + Qi ■ cf>i) + F(xi, 0J] 



(0) 



where: 



F(x l ,cj, t )=Y / \og{l±exp[~i( P]k -x l + q r ct> l )]} ±1 (7) 

jk 

In the above and following equations the upper sign ap- 
plies to fermions and the lower to bosons. By taking the 
continuum limit of the sum over phase space cell: 



(S) 



inipj-Xi+cij-tp^ 

(9) 



The density of states fli can now be used to obtain 
physical quantities of interest, for instance the primary 
(i.e. directly emitted by the hadronising source and not 
by subsequent hadronic decays) average multiplicity of the 
j th hadron species. Since every multi-hadronic state in the 
cluster has the same probability, it can be easily proved 
that this can be derived multiplying by a fictitious fugacity 
Xj the exponential factor exp[— i (pjk • Xi + • <pi)] for all 
k in Eq. (||) , taking the derivative of log fli with respect 
to Xj and finally setting Xj = 1: 



log fli (Xj 



(10) 



Most often, the physical quantities to be compared 
with experimental measurements are not those relevant 
to individual clusters, rather global ones, that is summed 
over all clusters in the event. This implies that sums over 
clusters with different four-momenta, quantum numbers 
and proper volumes must be performed and we are then 
led to consider all possible cluster configurations in terms 
of P, V* and Q and their probabilities. In the most general 
picture, one envisages the formation of a variable number 
N of clusters with probability Pn- For a fixed N, there 
shall be a conditional probability /*({Pi, Qi, V*}) rele- 
vant to the configuration {Pi,Qi,V7} of four-momenta 
Pi, quantum numbers Qi and proper volumes V* (i = 
1, . . . N); of course, any configuration must fulfill conser- 
vation laws, i.e. J2 i Pi = P and J2i Qi = Q where P and 
Q are the initial four-momentum and quantum number 
vector respectively. The corresponding distribution func- 
tion in the variables {P,, Qi, V{\ is: 



pO N pO 

/({P, Q ( , Vt}) = /*({P, Qi, Vi-j-}) J] -L (11) 



where M, 



ip} is the mass of the i cluster and P„° its 



energy. Thus, the overall average primary multiplicity of 
hadron j should be written as the sum of average primary 
multiplicities of single clusters (see Eq. (|l(J)), in a fixed 
configuration, weighted by the configuration probability 
/ and summed up over all configurations: 



N . 

= e p * [ n / d4p * dv * e ] f^ p > q 

N z=l Q, 
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where the symbol [Qi=i J d 4 PidVi SqJ stands for a iV th - 
uple integration and sum. The function / is in general 
unknown and depends on the dynamics of the cluster 
formation process. It can be further expanded according 
to the well known conditional probability decomposition 
P(AB) = P{A\B)P{B): 

f({Pi, Q„ Vi}) = <?({P, Qi}\{Vi}) H({Vi}) (13) 

where H({Vi}) is the probability distribution for the vol- 
umes {Vi} and g({Pi, Qi}\{Vi}) the conditional probabil- 
ity distribution of four-momenta and quantum numbers 
{Pi, Qi} once {Vi] are fixed. A considerable simplification 
would occur if g were equal to the w function below: 

w({P i ,Q i }\{V i }) P = 



6 4 (P-EiPi) 6 Q>SiQi l] 


^(jf)A(Pi,Qi,^) 


'u k E Q J^p k e(p°)n k 


5\P - S k P k ) 5 Q , SkQk 



(14) 

where 9 is the Heavyside step function. Because of the 
identity: 



x exp [i (P • x + Q • 0) + F(:r, 0)] (18) 

where: 
F(a:,0) = 

= |J-^V(2J i + l) / d 3 p logllie-'fe'^'^]^ 
{2ttY ^-r 1 J 

3 

(19) 

The fact that Q is the density of states of the EGC can 
be proved in an elegant way by showing that the w's in 
Eq. ( pT[ ) are just the probabilities of getting a set of four- 
momenta {Pi} and of quantum numbers {Qi} if a cluster 
of volume V, four-momentum P and quantum numbers Q 
is randomly split into N sub-clusters with fixed volumes 
{Vi} such that V = $^V<, for they maximise total en- 
tropy. We stress that the additivity of volumes in Eq. ( pj| ) 
applies to V^'s and not to proper ones V*'s because the 
splitting ideally takes place with no spacial overlap be- 
tween sub-clusters and Lorentz contraction must be taken 
into account. Let now w({Pi, Qi}\{V.}) p be such unknown 
probabilities and let us calculate the probability p of a full 
microscopic state assuming, for sake of simplicity, that 
{Pi} are discrete variables. According to the basic law of 
statistical mechanics, p turns out to be: 



i) 



N 

i— 1 "* 



N 



— io g n^(Ai) 



A,=l 



A j = l 



(15) 



the substitution of g in Eq. ( |T3| ) with w in Eq. (|T4|), turns 
Eq. @ into: 



(nj) =E P ^[I1 / dy^^D^log^A,) 

N i—l •* 

where Q{\j) is defined as: 



(16) 



= [IIE / ^P^P?)^)]^-^) S QtSiQi (17) 

i=l Q, 

It is a remarkable fact that ft(l) = Q is exactly the density 
of states of a single large cluster, here defined as equivalent 
global cluster (EGC) with four-momentum P = X^Pj, 
quantum numbers Q = Qi and volume (in the refer- 
ence frame where four-momentum is P) V — Vi and 
can thus be written in the same fashion as Qi in Eq. (^): 



p-^w^n ^t^ (20) 

where w({Pi, Qi}\{Vi}) P vanishes if P ^ P or Ff < 
or Q 7^ Si Qi and: 



^ w({P i ,Q i }\{V i }) P = l 
{PuQi} 



(21) 



The entropy S = — ^plogp should be calculated by sum- 
ming over all possible configurations {Pi, Qi} fulfilling the 
constraints on total four-momentum and quantum num- 
bers, and, once a configuration is fixed, over all possible 
microscopic states of the clusters. Therefore: 



N 

s = E w U E 

{Pi,Qi} i=l states. 



a 4 (P - /',>;^q q 



x log 



" 5 4 (Pi - P i;s )<5 Qi)Qi;s 



— 53 w 
{A,Qd L 



JX-(5 4 (P-P. s )Jq.q. 
logw + logj]^ 1 ' QQ 



(=1 



(22) 



where the arguments of w and Qi are implied. In the above 
equation advantage has been taken of the fact that log- 
arithm's argument is actually independent of the micro- 
scopic states of the clusters. In order to determine the w's, 
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S must be maximised with respect to all of them with the 
constraint = 1. This can be done by means of the 

Lagrange multiplier method which leads to the equation: 



as 



N 



dw({Pi,Qi}\{Vi})p " 

JV 

= -1 - logw - log Y[S 4 (Pi - P :s ) + log JJi7i + n = 

i=l i=l 

(23) 



implying that: 



w oc 



(24) 



The above equation, after a due normalization and tak- 
ing into account that the w's must vanish if P ^ T"V p 
or PP < or Q ^ J2i Qii coincides with Eq. (jig) , and 
this proves our statement. A different proof based on a 
direct calculation starting from Eq. ( |l7| ) can be found in 
Appendix A. 

By using Eq. @ with Aj = 1, Eq. @ can be written 

as: 



much fewer parameters, only those of the EGC (its total 
four-momentum, volume and quantum numbers). 

It is apparent from the previous derivation that the 
equivalence between a many-cluster system and one global 
cluster ultimately rests on the occurrence of configura- 
tion probabilities (P3). Furthermore, it can be realised 
by inspecting Eqs. Ql2^p"7|) that the equivalence holds in 
general for any observable A (not necessarily a Lorentz 
scalar), which can be written, for a given event configura- 
tion {P i; Q i: Vi}, as follows: 



A({P h Qi,Vi}) = 



c{i\Z 1 n t (p 1 ,Q t ,v l [ 



nili ni(Pi,Qi,Vi) 



(28) 



where £ is a linear operation, e.g. derivation or integra- 
tion. If this is the case, the observable averaged over all 
configurations with the probabilities (H) reads: 



(A) 



dV h(V) 



c(n(p,Q,v)) 



(29) 



Q{P,Q,V) 

Amongst such observables, one of the most important and 
general is the multi-hadronic probability distribution. For 
a given event configuration {Pj, Qi, Vi} it reads: 



w({Pi,Qi}\{Vi}) P = 

_ 5 4 (p - s^) s Q . SiQl Hi eyfMPi, Qi, SiVi 



n(P,Q,v) 



25) 



and Eq. (ha) as: 



f N f 

= J2Pn dV {[[ Wl\H({Vi})5(V - ZiVi) 



N 



x — \ogn(X jt P,Q,V) 



J2 p N f dVh N (V) ■^-logf2(X j ,P,Q,V) 

N J 



(26) 

A,=l 



where h N (V) = [JJ, J dVi\H({Vi})6(V - S l Y % ) is, by defi- 
nition, the probability distribution for the EGC with vol- 
ume V split into ./V sub-clusters. Since there is no further 
effective dependence left on N in the integrand of Eq. (Pq) , 
one can define h(V) = YIn •fW'iJv(V') so that the aver- 
age primary multiplicity of the j th hadron species finally 
reads: 



P(N 1 ,...,N K ) 
' 1 K 

— TT 



3=1 



dX-i 



N 

n 

i=l 



^({Aj},P«,Q«,yQ 

n l (p l ,Q l ,v i ) 



(30) 



In the above expression K is meant to be the total num- 
ber of hadron species; Xj are complex variables and the 
integration is taken on a closed path around the origin; 
the function f2i({Xj}) is a generalisation of the previously 
used f2i(Xj) with the insertion of fictitious fugacities of all 
hadron species at the same time. In fact, Eq. ([30|) has the 
form required in Eq. ( p8|) for the equivalence with EGC 
to apply. Eq. (^o|) can be obtained by inverting P's gen- 
erating function G(Ai, . . . , Xk)'- 



n u ...,n k 



E [E^ tate )| fix 



Ni ..... Nk stales 





■ ■ A K 


,N K )- 


Af 1 . 


\Nk 
■ A K 



= £p(state) Af . . . A^= J] m $p*£$r (31) 



states 2—1 

The last equality follows from: 



(27) 



dV h(V)— log f2(X J ,P,Q,V) 



This a noteworthy result because multiplicity in Eq. (|27J) 
no longer depends on the configuration of all clusters in 
the event nor on their number. Instead, it depends on 



K k n , 4 , p 

P( S tate)nAf =IK'II 

3=1 i=i 



P: 



5 Q,,Q !S 



i=l 



r?i(p i ,Q i ,F i ) 



(32) 



which is to be worked out thereafter as in Eqs. 

Besides the reduction in the number of parameters, 
the equivalence with EGC has another attracting feature: 
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for the canonical formalism to be a sufficiently accurate 
approximation, only the EGC has to be large, while there 
is no need to enforce a similar request for each individual 
cluster. Therefore, one would be allowed to treat canon- 
ically (as long as cluster-integrated quantities are con- 
cerned) even hadronising systems in wich single clusters 
are a priori known to be too small for a canonical treat- 
ment to apply individually. In fact, the transition from mi- 
crocanonical to canonical ensemble is based on an asymp- 
totic expansion of Q for large values of EGC volume and 
mass, through the saddle-point method, in which only the 
leading order is retained. To show this, first a rotation 
z — ix in the complex hyperplane for the four-dimensional 
integration in Eq. (fil]) is performed: 



lim 



-ioo+e i-tv in± 

/ Q 9 P-z+iQ-<M-F(-i 2: <M _ 



(27ri) 4 J_^ + ; J_ n (2tt)* 

= lim — — / d 4 z exp [P-z + log Z(z, Q)] (33) 
e^O (2tti) 4 J_ ioo+E 

where: 

Z(z, Q) = -i- f "d"0 exp [i Q • + F(-iz, <f>)} (34) 
(2tt)" J_„ 

One can recognize in Eq. (p4|), by looking at Eq. (|l^) for 
F, the expression of the canonical partition function flH, 
[l6| of an ideal hadron gas calculated for a complex four- 
temperature z. Now, the saddle-point asymptotic expan- 
sion of ft in Eq. ([33]) can be performed which, at the lead- 
ing order, reads: 



~ exp [P •/? + logZCS, Q)L 



(27r) 4 detH(/3,Q) 



(35) 



where (3 is a four- vector such that: 



d 



P ■ z + logZ(z,Q) 
P ' I+ 9^ l0gZ(z ' Q) 



z=0 



= 



M = 0,l,2,3 (36) 



and H is the hessian matrix 9 log Z/dz^dz 1 ' calculated for 
z = (3. Since F{— iz, <fi) is real valued (see Eq. (|l9|)) if z 
is a real four-vector and log Z(z, Q) is real too, according 
to Eq. d34|), then Eq. ( f36| ) states that /3 must be a real 
four-vector. Moreover, it has to be a timelike vector for 
the momentum integration in Eq. (]l9[) to converge. If /3 
is identified with the four-temperature vector (1/T)u, u 
being a unit timelike vector, the second term in Eq. ( |36"| ) 
is nothing but the expression, with negative sign, of the 
mean energy-momentum of a canonical system with four- 
temperature (3. In other words, the temperature can be 
defined for a microcanonical system (at the lowest order of 



an asymptotic saddle-point expansion) equating the given 
total energy to its expression for a canonical system having 
same volume V and the same set of quantum numbers Q. 
The canonical expression of entropy S — P ■ (3 + log Z can 
be recognized in Eq. ([35]) so that the well known Boltz- 
mann formula is recovered: 



Q oc exp S 



(37) 



The square root factor in Eq. ([35|) can be neglected in the 
derivation of most physical observables, such as multiplic- 
ities, because it depends on a fractional power of the par- 
tition function Z whilst the first factor is an exponential 
of it. Altogether, the expressions of multiplicities, multi- 
plicity distributions etc. in the canonical ensemble can be 
recovered. The canonical partition function ( |34| ) can also 
be written as follows: 

Z(J3,Q,V)= J2 e-^-^tW^ 

states 

= W d 4 Pe-^ 4 (P-P stat0 ) <5 Q , Qstatc 

states 

= J d 4 P c-P- p n(P, Q, V) (38) 

As has been mentioned, the volume V argument of fl in 
Eq. ( |38|) is measured in the frame where four-momentum 
is P. Indeed, there is a subtle difference between what is 
meant as proper volume in canonical and microcanonical 
ensemble. In the microcanonical ensemble the definition 
is clear because all states have a definite total momentum 
and the reference frame where it vanishes can be chosen. 
On the other hand, in the canonical ensemble, we can 
choose a reference frame where (3 — (1/T, 0) but this does 
not ensure that total momentum exactly vanish for all 
the states (it does so on the average with a small broad- 
ening around 0); hence, the actual proper volume (mea- 
sured with P = 0) does not coincide with the parameter 
volume used in the canonical partition function ( ^8|) with 
(3 = (1/T, 0). Nevertheless, if we think of a canonical en- 
semble as an approximation of a microcanonical ensemble 
in its rest frame with associated volume V* (according 
to Eq. (H^)), the four- vector (3 solution of Eq. j36| ) has 
vanishing spacial part. This means that the use of proper 
frame and proper volume V* in Q and Z in Eqs. (B5T) goes 
together with proper four-temperature (3 = (1/T,0). 

Two major points are worth being stressed. Firstly, 
as we have emphasised, the reduction to an EGC pos- 
sibly allows a canonical treatment for Lorentz invariant 
quantities even though physical cluster must be treated 
microcanonically and this implies that temperature may 
be a well defined quantity only in a global sense (at the 
level of EGC) while locally, at the level of single cluster, 
one should stick to microcanonically well defined concepts 
such as energy density. Secondly, a quantitative estimate 
of how large an EGC should be in terms of volume and 
mass for the canonical approximation to be satisfactory is 
highly desirable but it is not available by now. Neverthe- 
less, as the hadron gas has a huge number of degrees of 
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Fig. 1. Top: the typical configuration of the momenta of clus- 
ters in an actual high energy collision. Bottom: the configu- 
ration of the momenta of the same clusters originated from 
the splitting of one global cluster at rest. The two configura- 
tions are indeed equivalent for Lorentz invariant observables 
like hadron multiciplities. 



freedom, it can be reasonably expected that the validity 
of a canonical treatment should set in at relatively low 
values of volume and mass, though settling this issue def- 
initely requires very involved microcanonical calculations. 
At present, the legitimacy of the canonical approximation 
essentially relies on the agreement with the data. 



2.2 Back-boosting clusters 

As we have seen, the choice of probabilities (Q) for the 
configurations of the produced set of clusters is essential 
for the canonical approximation to apply. However, those 
configuration probabilities are unrealistic because hadrons 
emitted in a high energy e + e _ or hadronic collision should 
look like coming from one source at rest in the centre-of- 
mass frame whereas they typically emerge in a two (or 
more) jet- like structure (see Fig. |l|). Nevertheless, as long 
as one is interested in Lorentz invariant quantities such 
as particle multiplicities, it is possible to rearrange clus- 
ter momenta at leisure because there is no dependence 
on them as demonstrated by Eqs. (||,|l0|): the effective ar- 
guments of density of states and, consequently, average 
multiplicities, are mass and proper volume of the cluster. 
Then, the question arises whether it is possible to find 
a suitable rearrangements of cluster momenta so as to 
get expression for particle primary multiplicities equal to 
those obtained before in Eqs. ( 16 2q,|27|) for the EGC start- 



ing from a configuration probability distribution whatso- 
ever instead of the particular one (|l4|). 
This question may be restated more quantitatively as fol- 
lows: whether, given a general distribution of cluster four- 
momenta such as /* in Eq. (11), it is possible to find a 



" N 


fd 4 p 


_i=l 




" N 




_ i—l 


f d 4 P 



suitable rearrangement of cluster momenta in each event 
such that: 



f*({Pi,Qi,V*}) 0({Mi,Qi,V*}) = 
u4{Pi,Qi,V*}) 0({M u Qi,V*}) (39) 



where is a known distribution, M t = \fPl? is the mass 
of the i th cluster and O an arbitrary Lorentz invariant ob- 
servable dependent on cluster configurations. If the answer 
was affirmative, the task of reducing calculations of aver- 
age Lorentz invariant quantities to those of the EGC could 
be accomplished by choosing suitably related to con- 
ditional probabilities w in Eq. (14). Thereby, an effective 
cluster back-boosting from configurations as in Fig. [l] top 
to Fig. |l| bottom could be achieved. It must be pointed 
out, however, that /* and it* cannot be completely in- 
dependent of each other. In fact, since O in Eq. ( J39| ) is 
an arbitrary Lorentz invariant observable depending on 
{Pi}, the marginal distributions f*M({M i7 Q i; V*}) and 
u*m {{Mi, Qi, Vi}) obtained by integrating out the direc- 
tions of all Pi's in Minkowski space, must be equal, i.e.: 



/.Af ({Mi, Q,, V*}) = uMWi, Q*, V*}) (40) 

and this means that the actual cannot be completely 
arbitrary for the reduction to EGC to apply. This point 
will be discussed more in detail later on. 

In principle, the distributions /* and can be linked 
by a linear transformation through an unknown matrix X 
depending only on four-momenta directions P, P' in view 
of Eq. Mf): 



N 



i=l 



n ^ 



X(P i ,P i )u.({Pt 1 CU 1 V?}) (41) 



Provided that Eq. (^0|) is fulfilled, there are actually in- 
finitely many matrices X satisfying Eq. ( pd] ) and this en- 
ables us setting further requirements. For this purpose, it 
will be assumed that the dependence of X on Pj and P[ 
is realised only through Lorentz transformations L of the 
type (see ref. 0): 



R 3 (<p) R 2 (0) L 3 (0 



(42) 



transforming P[ into Pjj is a Lorentz boost along 

the z axis with positive hyperbolic sector £ and R/c are 
rotations around k th axis with angles ip £ [0, 2tt) and 9 £ 
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[0, 7r] 1^. Lorentz transformations of the type ( [42| ) allow to 
transform the time axis unit vector t — (1,0,0,0) into 
a given arbitrary timelike vector [fj"7| . Conversely, there is 
only one Lorentz transformation of this type transforming 
P[ into Pi (see Appendix B). Hence, Eq. ( pLT| ) is rewritten 
as: 



N 



i=l 



n / dL 



X({U}) ^({L^iP^Q^V*}) (43) 



where dL is the group measure for the Lorentz transfor- 
mations of the type (ff2|), namely [jlSj] : 



dL = 



d[L(i)ML(i)ML(i)} 



(44) 



Therefore, the problem stated at the beginning of the this 
subsection in Eq. ([39]) has been transformed into the quest 
of a solution X of the multiple-integral equation (ff3|), 
which is formally a Fredholm integral equation of the first 
kind. If a solution exists, then Eq. ( |39| ) is fulfilled for any 
Lorentz invariant observable O. It should be noted that if 
X satisfies Eq. @, then: 



n/dL 



X({U}) = 1 



(45) 



which can be obtained from Eq. fl43|) by enforcing the 
normalization constraints on both /* and it*: 



N 



i=l Qi 



({P i ,Q i ,V*}) = l (46) 



Henceforth, we will assume that a solution X of the Fred- 
holm integral equation ( p3[ ) exists and that it is a positive 
definite function. In fact, if X > 0, Eq. ([l3|) means that /* 
is a probability distribution of cluster four-momenta ob- 
tained by boosting clusters, primordially generated with 
a four- momenta probability distribution , with Lorentz 
transformations distributed with probability X. The re- 
quirement of positive definite X makes integral equation 
( 43|) not always solvable and it is not difficult to devise ac- 
tual examples. For instance, if is the distribution associ- 
ated with the splitting of an EGC into N sub-clusters and 
/* is a distribution characterized by constantly vanishing 
momenta of the clusters, no probability distribution X is 
able to transform u* into /* through (M3) because random 
boosts to the (small) non-vanishing momenta involved in 

yield non-vanishing momenta as well. 

In order to accomplish the reduction of calculations 
of average Lorentz invariant quantities to the EGC, we 



The axis index k here is not to be confused with the index 
i in following equations referring to the cluster which the boost 
is applied to. 



now have to specify the relationship between it* and w in 
Eq. wM- As a first step, let us define the new variables: 



M, 



(Pi) 1 P? 



M' 

V*4^ (47) 



along with the distribution u corresponding to in the 
V{ volumes (in analogy with Eq. (O)): 



PP 



^{PL^VlV^u^Pi^VlPf/M'^Xl^ (48) 



Evidently, V- are the cluster volumes measured in the 
frame where their four-momentum are P[ and the very 
reason of their introduction resides in their additivity, as 
it will become clear later. The distribution u can be ex- 
panded by means of conditional probabilities, in the same 
fashion as Eq. (|l3|): 



u({Pi,Qi,Vn) = v({Pl,Qi}\{Vl})Y({V;}) (49) 

In the second step, let us calculate the mean value of a 
general observable A by using Eq. (|4g|): 



(A) =J2 P " 



N 



fl^j^Pdv; 



i=l Qi 



f*({Pi,Q<,v*}) 



:A({P,Q i ,V*M l /P?}) 



E^ 



N 



N 

HE/ d^d^dln 



i=l Q 



X({U}) 



.({LT^P,), Qi, V*}) A{{P h Qi,V*Mi/P°}) (50) 



where P/v is the probability for N clusters to be produced 
in an event, (see Sect. 1). The variables defined in Eq. ( |47| ) 
and Eqs. ( (48| , [49| ) are now used to turn previous equation 
into: 



N 



N 

nE/ d 4 p'd^du 



i=l Qi 



x({u}) Y ({v;}) 



xv({Pi QMVI}) A({U(P/), Q h K'PfVU(P') }) (51) 

No weighting factor is implied in the P[ = L^ 1 (P^) trans- 
formation as | det Li \ = 1. If A is a Lorentz invariant ob- 
servable like O in Eq. (^), then it must be independent of 
Lj and the integration of X can be performed separately 
yielding 1 in accordance with Eq. (Eq). Hence: 



(o) =E P ^ 



N 



nE/ d4 w 



i=l Q 



v{{PiQl,}\{Vl}) 



xY({Vn)0({Ml,Qi,V?P!°/Ml}) 



(52) 
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The form of Eq. (^2|) suggests that one should set v = w. 
Indeed, if O can be written as in Eq. (p8|), the whole 
derivation in Subsect. 2.1 for average primary multiplici- 
ties can be repeated, leading to: 



(O) = Y,Pn 



N 



n£/d 4 p;dK 



Li=l Q 



'(.{Pt,Qi}\{V[})p> 



xY({V!})O({Ml,Qi,V!Pl /Ml}) 



Y({V{}) 0{M', Q, P'°SiVf/M') 



n / < 

<=i 

J dV y(V) 0(M', Q, V'P'°/M') 



(53) 



where O in the last expression is meant to be the same 
observable calculated in an EGC with four-momentum P' , 
mass M' = VP' 2 , quantum numbers Q = J^. Qi, volume 
V' = ^7 ( m the frame where four-momentum is P') 
and: 



AT 



N ■ 

n / ^ 



y(TOW--W) (54) 



Now we can finally write the sought relationship between 
and w on the basis of Eqs. ( |48| , [f9| ) and the condition 



«»({Pi,Qi,^*}) 



JV 



= w ( {Pi ,Q i }\ { V*^}) P ,Y({V*^})f[^ (55) 



Therefore, for any given /* configuration probability dis- 
tribution, the integral equation ( f43| ) for the reduction to 
the EGC finally becomes: 



f*({Pi,Qi,V*}) = 



n f**-^- 



*({U}) 



x«({L7 i (fl) 1 (^}|{Vrq 7 ^ 5 })/-nWq^ 



(56) 

with X function to be determined. Whilst w is given by 
Eq. (Q, Y is determined by integrating both sides of 
above equation in the variables P; 's and summing over the 
Qi's. After the change of variables Pi — > P[ in Eq. (^7|) , 
one obtains: 



whose left hand side is apparently the actual marginal 
probability distribution for cluster proper volumes. Thus, 
the function Y should satisfy the integral equation: 



H*({V*}) 



dP 



dV{ ...dVj 



N 



N 

II £ 

i=l Qi 



d 4 P' 



p/O 



xw ({p!,Q i} \{V*^}) P ,Y({V*^})(58) 



Like for Eq. (|43[) , we will simply assume the existence of 
a positive definite solution for Y. 

We are now in a position to answer the question 
whether one can make a reduction of Lorentz invariants 
calculation to the EGC starting from a /* distribution 
whatsoever by solving the two integral equations (|56"|,|58|) . 
According to the discussion with respect to Eq. (|4(f) , this 
is not true in general because the first integral transforma- 
tion through X does not affect the distribution of cluster 
masses. Likewise, Eq. ( p8| ) only implements the accordance 
between the actual proper volumes distribution and that 
determined by the required insertion of w, hence it has no 
effect on masses. Therefore, there must be a definite con- 
straint to be fulfilled by /* for the reduction to the EGC 
to be possible and this is concerned with the conditional 
probability distribution of cluster masses and quantum 
numbers once their volumes are fixed. This is best seen by 
integrating both sides of Eq (|5q ) with the usual change 
of variables Pi — > P[ in Eq. (|47|) and taking into account 
the normalization condition in Eq. (45): 



II Jd 4 Pis(Mi 



f*({P h Q h V*}) 



fMWu Q 4 , V*}) = gMWu Qi}\{V*}) H,{{V*}) 
f[f^P^S{M,-JPP 



xw ({ P !,Q i} \{V*^}) P ,Y({V*^}) 



(59) 



where /* has been decomposed into the product of a con- 
ditional probability distribution of masses and quantum 
numbers g*M for fixed proper volumes and the probabil- 
ity distribution of proper volumes P* . The above equation 
can be solved to determine g*M by using the expression 
of P* in Eq. (jH^) . Indeed, the insertion of P* does but to 
set the due overall normalization of g*M to 1. Thus: 



Li=l 

N 

|_i=l Q, J i 



f*({Pi,Qi,vn) = 



Mi 



™({Pl,Qi}\{V*^>})p'Y({V*-^}) 

i 

(57) 



9*M({Mi,Qi}\{V*}) 



N f 



■ Mi 



N 



xw({P!,Qi}\{V*^}) P < Y({V*^}) 



(60) 



where Nf is the normalization factor obtained by inte- 
grating the numerator in the ilp's and summing over Qi. 
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Since Y can be determined as a function of w and 
through the integral equation (58), EQ. (|6(]) states that 
g*Mi the actual probability distribution of cluster masses 
and quantum numbers for fixed proper volumes, ought to 
fulfill an independent constraint. This is definitely not a 
trivial requirement and must be taken as a further as- 
sumption to make it possible the existence of an EGC. 
Most likely, the actual g*u shall not be equal to the right 
hand side of Eq. ( |60|) but, hopefully, close enough to it to 
make the equivalence with the EGC a good approxima- 
tion, thus justifying the surprising accuracy of multiplicity 
fits performed under that assumption. 

It is of utmost importance to stress that the four- 
vector P' in Eq. (p3|), whose appearance is related to 
the particular form of w (see Eq. (fl4|)), is not the initial 
four-momentum of the collision. Instead, P' is the four- 
momentum of the EGC to be set as a result of cluster 
back-boosting and can be taken as a free parameter. For 
any chosen value of P', there shall possibly be two corre- 
sponding solutions Xpi and Yp> of the integral equations 
(^.^). There are, however, some restrictions: M' — V P a 
must certainly be smaller than centre-of-mass energy y/s 
and greater than the actual maximum value of J^. Mi, the 
sum of cluster masses. The latter constraint is the most 
relevant, for it might spoil at once our EGC reduction 
scheme falsifying the relationship (|60|). For instance, a g*M 
distribution can be envisaged which is in agreement with 
that needed for an EGC with mass M <C y/s except for a 
small tail or bump over which the sum of cluster masses 
is very close to y/s. In this case, it is apparent that the 
equivalence with an EGC having one definite mass turns 
out to be impossible. Nevertheless, in principle, it should 
be possible to split collision events into suitable subsets 
having different values for max[^ i Mj] and link them to 
different EGC's four-momenta P', just because P' is a free 
parameter. Pushing along this idea, one is led to consider 
a continuous distribution T(P') to be chosen so as to re- 
store the agreement between the actual and EGC-splitting 
distributions of J^. M* = M t : 



G(M t ) 
xS(M t 

= J d 4 p' r(p') 



N 

IIE 

i=l Q, 

SiMi) 

N 

IIE 

Li=l Qi 



d 4 P dV* 



f*({Pi,Qi,V*}) 



d 4 p'dv; 



S(M t - SiMl) 



xwdPiQMVlVp'YW}) 



(61) 



where Eq. (|5q) h as been integrated by using the change 
of variables ( |17| ) and the normalization condition (Eq). 
Again, we will assume that a positive definite solutions 
of this integral equation exists. Of course, the introduction 
of the function T(P') requires Eqs. ( |56] , |58| ) to be changed 
accordingly: 



f*({Pi,Qi,V*}) = 



N 



n dL 



Mi 



X({U}) 



J d 4 p> r(p')«;({Lr 1 (P0.Q<}IWL?^ 



Mi 



(62) 



and: 



H.{{V?}) = 



N 1 

IIE/ d4 ^ /d 4 P'T(P') 



Li=l Q 

Mi 



Mi 



x w ({P/,Q i }|{^-i}) P ,y({^- 7 l}) (63) 

i i 

Altogether, the equations ( |6l| , |6^ , |63| ) are a system of cou- 
pled integral equations for the unknown functions X, Y 
and T to be solved at the same time in order to achieve 
the reduction to a superposition of EGC's with different 
four-momentum and volume. Note that in Eqs. (pHp2Lp3|) 
X and Y have been assumed to be independent of P\ yet 
this dependence can be introduced without affecting most 
of previous and following arguments. 
As has been discussed with regard to Eq. (|o|), these equa- 
tions are not sufficient to ensure that any configurational 
distribution /* is equivalent to such a superposition, as 
far as the calculation of Lorentz invariant quantities is 
concerned. In fact, after the introduction of T, the con- 
straint (|o|) is left almost unchanged: 



9*M({M h Q i }\{V*}) = — 



f M' 

II ^n jskmi-M'i) 

/M' M' 
d 4 P' r(P')u,({P i ',Q i }|{TA*-^})p,y({^-^})(64) 

and should be fulfilled independently once Y and Y have 
been determined as a function of the given /* and w. Fi- 
nally, the average value of a Lorentz invariant observable 
O reads (see Eqs. d||E3)): 



(O) =E P ^ 



A? 



N ■ 

n / ^ 

i—l 



d 4 P' Y(P')Y({V?}) 



xO(M',Q,EiV^P' /M') 
= J dV J d 4 P' Y(P') y(V) 0(M', Q, V'P'°/M') 

/t M' M' 

dV* J d 4 P> Y(P')y(V*j^)0(M',Q,V*) 

= J dV* J dM' X (M',V*)0(M',Q,V*) (65) 

where V* is the EGC's proper volume and: 



X (M',V*)= / d 4 P'<5 



M' 



M' 



(66) 
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If the distribution x allows only large values of volume 
and mass, the canonical approximation of the last expres- 
sion in Eq. ( |65|) could be a satisfactory one. Specifically, 
the microcanonical average 0(M' , Q, V*) is to be replaced 
by the canonical average Ot(T(M', V*), Q, V*) calculated 
for the suitable temperature T dependent on mass and vol- 
ume (see discussion at the end of Subsect. 2.1), so that: 



(O) = J dV* J dM' x{M', V*) 0(M', Q, V*) 

~ J dV*J dM' X (M',V*)0 T (T(M',V*),Q,V*) 

f f dM 1 
= J dV*J dT X (M'(T, V*), V*) — O t (T, Q, V*) 

= J dV*J dT C(T,V*)O t (T,Q,V*) (67) 

Thereby, the actual average (O) has been written as a 
superposition of canonical averages at different tempera- 
tures and volumes. The last integral in Eq. ( |67| ) can be 
replaced with the mean value of the observable calculated 
for average values of T (= T) and V* (= V ), namely: 

J dV* J dT((T, V*)O t (T,Q,V*) = O t (T,Q,V*) 

(68) 

In principle, T and V depend on the particular observ- 
able O (e.g. on the hadron mass if O is a multiplicity). 
Notwithstanding, in all previous studies g,f§| on aver- 
age multiplicities of hadrons in the canonical ensemble, 
the single temperature and volume average in the right 
hand side of Eq. (^) has been used, tacitly assuming that 
T and V were independent of O. This is a reasonable 
assumption only if £ is a strongly peaked function of T 
and V* or, in other words, T and V have small fluctua- 
tions. On the other hand, one of the main results of such 
analyses is indeed the very good agreement between the 
data and single-temperature and volume fits at a fairly 
constant temperature value, which somehow a posteriori 
justifies the above assumption. That finding also provides 
a strong indication that hadronisation occurs at a critical 
value of cluster energy density g{ . 

Instead of deriving canonical formulae as a limiting 
case of microcanonical ones, it might be possible, in prin- 
ciple, to enforce the exact validity of the canonical en- 
semble by a suitable choice of T, that is of EGC's mass 
fluctuations. For this purpose, T must be set equal to the 
probability density cr(P') of getting a four- momentum P' 
in a canonical system with four-temperature j3, volume 
V' = "Yin V( and quantum numbers Q: 



0-(J y )/J,Q,V = 



Estates 5\P'- Pstate) e^ P ^S Q . Qa 



Z(f3,Q,V) 



where Z is the canonical partition function. Note that in 
Eq. (69) the EGC's four-momentum distribution function 
depends on V — J2 { V( besides P', but this does not affect 
any of the previous and forthcoming arguments provided 
that integrations are performed in the appropriate order. 
An integral equation corresponding to equation (|6l]) in the 
canonical ensemble can be obtained by replacing T(P') 
with a superposition of <j(P') ! s for different temperatures, 
with = (1/T,0): 



T{P')= / dT <^(T) cr(P')p,Q.v 
Jo 

= n(P , ,Q,V') I"" dT 1 
Jo 



Z(T,Q,V) 



-P'°/T 



(70) 



where W{T) is an unknown temperature distribution. By 
plugging Eq. ( |70| ) into Eq. (|(3l]) one obtains again a Fred- 
holm integral equation of the first kind for <?": 



JV 



5{M t -SiMi) 



G(M t )=[ dT <P(T) f[J2[ d 4 P/dT// 
Jo Li=l Qi J 

' nwi) (7i) 

which, however, may not always have a solution for the ac- 
tual G(M t ) might vanish if M t is less than a finite thresh- 
old value whereas the right hand side does not, regardless 
of t^/'s shape. This is indeed one of the key properties of 
the canonical ensemble: for any temperature, the distri- 
bution of EGC's invariant mass (and of SMi too) always 
has non-vanishing tails down to zero and up to infinity 
whilst a superposition of EGCs with definite masses may 
not. Therefore, the sought exact equivalence between the 
actual system and a canonical superposition of global clus- 
ters with finite temperature and volume has one more dif- 
ficulty with respect to the microcanonical superposition. 
It is interesting to derive the average value of the Lorentz 
invariant observable O in this case similarly to what has 
been done before in the microcanonical case. This can be 
done by replacing T(M) in Eq. © with Eq. @ : 



Y({Vn) 



p'O 

OiM'^—SiV!) 



(69) 



pOO JV p 

(0)=Y,Pn / dT«^(T) J[ dV> 

x /d 4 P'c- p '° /T Q ( P '^^ V D 
J Z(T,Q,£iV?) 

= ^°dT^(T) JdV'y(V) Jd 4 P' c- p ' a ' T 
x ^'^ 0(M',Q,V'P'°/M') 
= J dT<P(T) J dV y(V')0 T (T,Q,V') (72) 
where y(V) is the same function defined in Eq. (|54| ) and: 
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T (T,Q,V') = 

]lo ,e- pm / T a(P',Q,V) 



d 4 P' 



Z{T, Q, V) 



0(M',Q,V'P'°/M')(73) 



is apparently the average value of the observable O for an 
equivalent global cluster at temperature T, volume V and 
quantum numbers Q. 



2.3 Summary 

Before moving to the central topic of this paper, i.e. tran- 
verse momentum spectra, it is worth summarizing the 
main points of this Section. 

— The statistical hadronisation model assumes hadrons 
to be produced from a set of clusters with dynamically 
generated configurations in terms of quantum num- 
bers, volumes, mass and momenta; within each cluster, 
all hadronic states are equally likely. 

— Since cluster configurations are dynamically generated, 
the statistical ansatz needs a supplementary dynam- 
ical information. However, Lorentz invariant observ- 
ables, such as hadron multiplicities, are unaffected by 
a change of cluster momenta. Then, in order to reduce 
the number of free parameters to a minimum, cluster 
momenta distribution can be modified so as to achieve 
an equivalence with the calculation for one global clus- 
ter (EGC); this is done through a system of integral 
equations. 

— Provided that positive definite solutions of the afore- 
mentioned integral equations exist, the reduction to 
the EGC ultimately relies on a particular form (see 
Eqs. ( |60| , |64| )) of cluster mass and charges fluctuations 
at fixed volumes. If this is not the case, the equivalence 
is spoiled. Still, the actual distribution could be close 
enough to that particular form to ensure a sufficiently 
accurate approximate equivalence. 

— If mass and volume of the EGC are large, the basic 
microcanonical treatment can be replaced by a more 
manageable canonical one. Thanks to the equivalence 
with the EGC, this should be possible even though in- 
dividual physical clusters are too small to be treated 
canonically. Therefore, temperature could be a well de- 
fined quantity only in a global sense and not at locally 
for each single cluster. 



3 Transverse momentum spectra 

We can now proceed the main subject of this paper, that is 
transverse momentum spectra. Unlike average multiplici- 
ties, this observable is not a Lorentz invariant one and a 
special treatment is necessary. To start with, let us con- 
sider the spectrum of the j hadron species relevant to 
the i th cluster: 



dpT 



Or) 



Pi,Qi,V" 



(74) 



The possible dependences of this spectrum on the four- 
momenta Pi are crucial for forthcoming arguments. If a 
Lorentz transformation A is performed and the new trans- 
verse momentum is p' T = A(p)y, the corresponding spec- 
trum must have the same functional dependence on p' T 
as the transverse momentum spectrum associated with 
the cluster whose four-momentum is Lorentz-transformed 
with the same A. Hence (Q, and V* have been omitted): 



Apz 



dA- 1 (p) T 



dpi 



dn 



j i 



dpi 



(pt) 



A(A) 



(75) 



A rotation around z axis does not change pr', thus, ac- 
cording to the above equation, the spectrum must depend 
on Pf and Pf only through P^ and not on the azimuthal 
angle. Furthermore, px is not changed by a boost along z 
axis either, so pji must depend on Pf and Pf only through 
the combination if 2 - Pf 2 = Mf + Pf T . Therefore, the 
spectrum ought to have a dependence on Mi and PiT only, 
so that: 



dpi 



{pt) 



Pt,Qi,V* 



dry 
dpi 



■(pr) 



(76) 



Indeed, what we are really interested in is the overall av- 
erage transverse momentum spectrum of a given hadron 
species (labelled by j). To obtain it, the distribution of 
mass, volume and quantum numbers of the clusters must 
be folded with the sum of the single cluster spectra for a 
given configuration, yielding (see Eq. j5C|)): 



{ dPT )~Y. p > 



N 



N 



i=l Q 



N 



E^Cpx) 



Mi,P iT ,Qi,V* 



/*({Pi,Qi,VT}) 



(77) 



The distribution /» can be replaced with the right hand 
side of Eq. (|6^) and, by using the change of variables (f|7j), 
Eq. fin\ ) becomes: 



/ An 3 \ = ^ 

N 



\dp T ' 



p. 



N 



N 



i=l Q 



mu}) 



J d 4 F' T(P') w{{Pi Qli}\{V!}) P ,Y{{Vl}) 



N 



dn, 



r-f dp T 



m,,u(p[)t,q,,v;^- 



(78) 



The physical meaning of this expression is essentially that 
the average transverse spectrum of the hadron j is the 
convolution of a distribution X of Lorentz transformation 
on clusters (assumed to be positive definite, see Sect. 2) 
with the spectra obtained from EGC splitting into iV 
sub-clusters. The integration over Li can be considerably 
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simplified by using an alternative decomposition of the 
Lorentz transformation L defined in Eq. ( |42] ) (see Ap- 
pendix B): 



Ufa) R 3 (<« ii(C) 



(79) 



where and R& are to be understood as in Eq. (|42|), <j) € 
[0, 2ir) is an angle and r] £ (— oo, +oo) and ( £ [0, oo) are 
hyperbolic sectors. This particular decomposition, along 
with its associated measure (see again Appendix B): 



dL = — sinh 2£ drj dip dC 



(80) 



is advantegeous in that it allows to integrate away at once 
two parameters in Eq. (|7"8|). In fact: 

U{Pl) T = (L 3 (7 ?i )(R 3 (^)(L 1 (Ci)(^))))T = 

= (R 3 (0 4 )(Li(G)(^)))t = UCm)T (81) 

because both L 3 and R 3 leave the transverse component of 
their argument unchanged. Thereby, the transverse com- 
ponent of the four-momentum P[ undergoing a Lorentz 
transformation like (^2|) or (|7^), is actually the same ob- 
tained by applying one suitable Lorentz boost along x 
axis. The parameter Q > is related to the four-velocity 
associated with this transverse boost, namely: 



sinhCi = flTilTi = u T i cosh(i = = 



while the associated measure reads: 



- sinh 2£i dC, = u Tl du Ti 



Therefore, Eq. (|7q) becomes: 



/dnj_ 
\dpx 



E^ 



nE/ d4 w 



du 2 Tl 



i=l Q 



1 + u 2 Ti 

(82) 

(83) 
Xt{{u 2 Ti }) 



J d 4 P' n^wmMmvp'YW}) 

o (84) 



N 



Edriji 
d^ {pT) 



i=l 
,2 



M i M(uT i )(P!)T,Qi,V'-£- 



where Xxdu 2 ^^) is defined as: 



X T ({u 2 Ti }) = 



n / a ^ / d 

i=l J-oo JO 



X({U{rii,(j>i,UTi)}) 
(85) 

Before doing anything else, it is advantageous to work out 
the integral: 



in which the other arguments of driji / dpr have been omit- 
ted. In fact, at this stage, we can take advantage of a 
physical information concerning the range over which the 
distribution Xt is expected to be significantly different 
from zero. Indeed, it is well known that in high energy col- 
lisions hadrons are mostly emitted with a limited trans- 
verse momentum (with respect to event axis in e + e _ or 
beam line in hadronic collisions) of the order of few hun- 
dreds MeV, slowly increasing with centre-of-mass energy. 
Consequently, in the framework of a cluster hadronisation 
model, most clusters must have a little transverse momen- 
tum and, by taking as educated guess their mass values 
between 1 and 3 GeV, transverse four- velocities generally 
< 1. In other words, the distribution Xt(u 2 ^ 1 , . . . , iij~ N ) is 
expected to be significantly different from zero only in the 
region v£ ri < 1 i — 1, . . . , N over a reasonably large centre- 
of-mass energy range. Certainly, at very high energy, the 
increase of transverse phase space (for instance, the in- 
crease of radiated gluon pr in e + e~ collisions) is likely to 
extend the effective transverse four-velocity domain be- 
yond 1. Nevertheless, in the present work, we will confine 
our attention to a centre-of-mass energy range where the 
condition <C 1 applies. Under this circumstance, one is 
allowed to expand all driji /dp^s in powers of r a\ i starting 
from a suitable point u 2 ^ common to all of the clusters: 



-(pt) 



dpr 



Li(« Ti )(P/)t 



T ji( u Ti) — T ji( u T) 



u 2 T ) + 0((u 2 Ti - u 2 ,^ 2 



(87) 



so that the integral ySq ) becomes: 



N 

n 

i=l 



°du 



Ti 







X T ({u T i}) 



N 



E r ^(4) 



Bt- 

1 2 ~ { 



+j£r (4) (4* - 4) + o{{u 2 Tl - u 2 T f) 



The Xt distribution will be assumed to be such that even 
the first order term in the expansion ( E37j ) can be neglected. 
Then, only the zeroth order term will be retained and 
Eq. ( p4|) turns into: 



JV 



d 4 p' r(p') 



xw({Pl,Qi}\{Vl}) P >Y({vn) 
-(Pt) 



N A 



i=l 



dp 7 



M I X 1 (G T )(P/) r ,Q,,V;' li J- 



(89) 



N 

n 

i=l 



°°du 2 Ti 



X T {{u 2 Tl }) 



N A 



dpi 



(Pt) 



where ut = w u, T . 
Now the reduction to the EGC may apply due to the fact 
Li(u Ti )(P!) T that all clusters in the above equation are boosted by the 
(86) same transverse four-velocity ut along x axis. In other 
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words, the sum of transverse momentum spectra of all 
clusters for a given configuration of four- momenta and vol- 
umes, with subsequent integration over them weighted by 
w, must yield the transverse momentum spectrum of the 
given hadron in an EGC boosted with four- velocity ut- 
This statement can be proved not only for transverse mo- 
mentum spectrum but, more generally, for the invariant 
four-momentum spectrum. For this purpose, the convolu- 
tion equation Eq. ( |S9l) is generalised to four-momentum 
spectra: 



The proof is similar to that for average primary multi- 
plicities and can be conveniently carried out by using the 
expression of fli in (|l^) with the obviously needed re- 
placement, for single-mass-valued particles (i.e. not reso- 



d 3 p = / d 4 p S(p° - a/p 2 + to 2 ) 



(95) 



A' 



i=l Qi 



xv>({Pt,Qi}\{Vft)p'Y({Vft) 
-(P) 



(90) 



To complete the proof, four-momentum spectrum for a 
given configuration must be expressed as in Eq. (p8|), ac- 
cording to what has been established in Sect. 1. Since: 



drio 



■(P) 



d 4 p 
so that Eq 



Li(« T )(F i '),Q i ,V i * d 4 p 

can be rewritten as 



P'.Qi,V* 

(91) 



\d 4 pl 



N 



N 



II E/ / ^p'np') 



xw{{PiQi}\{Vl})p>Y{{Vl}) 



N 



£^fr l (*r)(fO) 



i=l 



d 4 p 



P!,Qi,Vj- 



(92) 



Now, the contributions of primary hadrons, i.e. directly 
emitted from hadronising clusters, are separated from that 
of decay products of heavier hadrons: 



/ dflji \ _ I drijj \ Primary 
\d 4 p/\ d 4 p I 



E 



/dnj£\ 
\ d 4 P I 



(93) 



where dn,ji/d p\ i is the four-momentum density of 
hadrons j stemming from the decay of the primary hadron 
k either directly or through intermediate steps (i.e. k — > 
h — > ... — ► j). The primary part of the spectrum can be 
calculated by using a technique similar to that for average 
primary multiplicities, introduced in Eq. ([To|). Instead of 
taking the derivative of log Qi with respect to a fictitious 
fugacity Xj, one has to take the functional derivative of 
log fii with respect to a fictitious 'fugacity' function Xj (p) 
depending on four-momentum p and set to 1 thereafter: 



dflji primary 5 log Q l [Xj] 



d 4 p 



5X 3 {p) 



(94) 



\j(p)=\ 



It is worth emphasizing that, unlike the derivative with 
respect to Xj, the functional derivative does not re- 
turn a Lorentz invariant quantity. The secondary spec- 
tra driji/d 4 p\ k ^ : > can be expressed as a convolution of 

the 

four-momentum density of hadron j stemming from the 
decays of primary hadrons k with four-momentum t) ei- 
ther directly or through intermediate steps, with the four- 
momentum density of primary hadrons k expressed again 
as a functional derivative: 



at i tra ariji/a p\ J can be expressed as a convoiutio 

JI / d 4 P/dV/ / d 4 P' T(P') a cluster-independent kernel function Kj k (p,t) (i.e. 
n J J four- momentum density of hadron j stemming from 



dn 



31 



d 4 p 



dHKjkip, t) 



8\ogQi[\ k ] 



5X k {t) 



(96) 



\ k (t)=i 



Hence, the four-momentum spectrum of the hadron j for 
a given configuration of clusters, i.e. the last factor on the 
right hand side of Eq. (B5), reads: 



A 



dn 



+ J2 ! dHK jk (L^{u T )(p),t) 



N 

E 



8 log Qi [Xj 



6 log Qi[X k ] 



5X k (t) 



SXj(l^(u T )(p)) , 

s 



£ t- 1 (« T )(p) 



E 



AT 



8X k (t) 



\i=l 



(97) 



where the final setting Xj — 1, = 1 are implied. The 
operator between square brackets in the above equation, 
defined as Ct-i,~ u \, is a linear one, so that the four- 

momentum spectrum of the hadron j for a given configu- 
ration of clusters has the form required in Eq. (^8|) for the 
reduction to the EGC to apply. Therefore, using Eq. (|97|) 
and proceeding like for Eq. (p3), the four-momentum spec- 
trum in Eq. tod) can also be written as: 



N 



N ■ 



Y({Vf}) / d 4 P> T(P') 



/2(P',Q,A^/) 
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dV d 4 P' r(P')y{V 



(/?(p',q,v")) 



^(P',Q,y) 



(98) 



From now on, calculations will be carried out in the 
canonical ensemble only. Formally, this amounts to take 
the T function as in Eq. (^0|); we have mentioned in Sect. 1 
that this choice might not be an appropriate one if a lower 
bound on the sum of cluster exists, yet we will neglect this 
possibility. Moreover, we will assume a temperature distri- 
bution function ^/ (see again Eq. (|70|)) equal to, or at least 
very close, to a Dirac's delta or, in other words, a unique 
temperature (see discussion at the end of Sect. 2). These 
assumptions amount to enforce an EGC's four-momentum 
distribution function Y in Eq. ( p8|) equal to a in Eq. ( |69|) 
with (3 — (1/T, 0). The main justification of it is a very 
strong indication in favour of a critical come out from the 
analysis of multiplicities. Whether the same critical value 
is retrieved in the analysis of transverse momentum spec- 
tra, it is just the main issue to be studied in this paper. 
Hence: 



[ dV> [ d*P> y(V>) ^Q.n e -P-/T 
J J yv ; Z(T,Q,V) 



dn 3 
d 4 P 



n(P',Q,v) 



dV y(V) 



Z(T,Q,V) 



-I 



dV y(V) 



Z(T, Q, V) 



(99) 



where advantage has been taken of the linearity of C and 
Eq. (HI) has been used. The integration over volumes V 
weighted by the y(V) distribution is now replaced by a 
mean volume V which is assumed to be independent of 
hadron species; this would be an exact operation were 
not for th e de pendence of chemical factors on the volume 
(see Eq. (104) below). Thus, the final result is the four- 
momentum invariant spectrum of hadron j in a canonical 
EGC with volume V and temperature T evaluated at the 
four- momentum Li 1 (ut)(p), that is the same spectrum, 
evaluated at the four-momentum p, in a canonical EGC 
boosted with Li(«t) (see Eq. (pl|)): 



\d 4 p/ ~ Z(T,Q,V) 

= C p (Z(^ fiT ,Q,F/7T)) 
Z(T,Q,V) 



where = L 1 (m t )((1/T, 0)) = (j T /T, j T p T /T, 0, 0) 



I + u 2 T /T,ur/T, 0,0). 

The transverse momentum spectrum can be obtained 
by integrating the above spectrum in mass, p z and tp in 



cylindrical coordinates. In fact, it can be shown that this 
operation, from a formal point of view, amounts to a redef- 
inition of the fu nctio nal derivative and integral operator 
in Eqs. @@|§,|0^): 



/dnj_ 
\ dpT 



k 



5logZ([\j],j3u T ,Q,V/jT) 



6\j{p T ) 



d 4 tK Tjk ( PT ,t) 



A 3 - = l 

51ogZ([A fc ],/3 UT ,Q,F/7T) 



sx k (t) 



(101) 

A fe =l 



where Kxjk is the transverse momentum density of hadron 
j stemming from the decays of hadron k with four- 
momentum t. We are now going to show how to perform 
a detailed calculation of the above expression for different 
temperatures and transverse velocities. 



4 Calculation of primary and secondary 
spectra in the canonical ensemble 



Before working out Eq. (|10l| ), it is worth summarizing the 
key assumptions that have led to that expression: 



1. 



the reducibihty of the actual set of clusters to an EGC 
as far as the calculation of mean values of Lorentz 
invariant observables is concerned; this has been ex- 
tensively discussed in Sect. 2 and requires that inte- 
gral equations (|5q, |58|) for fixed EGC mass, or integral 
equations ( |6l| , |62| , |63p or variable EGC mass, have pos- 
itive definite solutions; 

an EGC large enough to allow the use of canonical 
ensemble; 

transverse four-velocities ut = PtIt of the clusters 
much smaller than 1, making the replacement of a 
transverse four-velocity distribution with an average 
value ut a good approximation. 



Eq. (101) has two terms which are the contributions to 
the final spectrum of the primary and secondary (i.e. pro- 
duced by hadronic decays) hadrons respectively. Taking 
into account that the canonical partition function reads 



Z([{A fc }],/3 fiT ,Q,^) = e^+^IlM]^) 

IT (27T) 7 V_ 7r 

(102) 

where: 



(100) F c ([{Xj}],<Pi) = 



V 



7t(2tt)3 



E( 2J , + 1) 



x J d 3 p log [1 ± \{ PT ) exp(-/3 Sr ■ Pj + iqj ■ 0)] ±1 (1O3) 

the first functional derivative in the right hand side of 
Eq. (|j~0l|) yields: 
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, dnj > primary = F(2J 3 + 1) +1 



-OO p TV 



' ' Z(T, Q, V) 



OO J —TT 



(104) 



with 7r = y 1 + defined below Eq. (U00|). The ra- 

tios of partition functions in the above equations have 
been defined as chemical factors in refs. [UP as they con- 
tain the dependence of the production rate of the consid- 
ered hadron on its quantum numbers. The integration in 
p z and ip can be done analytically, yielding: 



/dnA P" mar y V(2Jj + 1) ^ 

(t^/ r = V(Tl) + "irpr 

Vd ^ 7 2^^1+1^ 



T 



V T / Z(T,Q,V) 



(105) 



p^, + m 2 is the transverse mass and Ki, 



where tut = 

Io are modified Bessel functions. By integrating the above 
spectrum, the formulae for primary multiplicities quoted 
in refs. can be recovered. At a constant tempera- 

ture of about 160 MeV found in the analysis of multiplic- 
ities, the effect of quantum statistics can be neglected for 
all hadrons e xcep t pions and only the first term of the 
series in Eq. ( |l04| ) can be retained. This approximation 
makes the spectrum shape unaffected by the chemical fac- 
tor Z(Q — qj)/Z(Q) which can be embodied in an overall 
normalization factor. Thus, in the Boltzmann limit: 



/ drij \ Primary 
\ dpT 



xKi 




Tot pt 



where 



(106) 



is the primary multiplicity. For resonances 



with wid th _T > 1 MeV, the single mass value mj in 
Eq. ( |105| ) is replaced with a relativistic Breit-Wigner dis- 
tribution over an interval [m — Am,m + Am], where 
Am = min[m - mthreshold, 2P] and TOthreshold is the mini- 
mal required mass to open all decay modes. 

The calculation of the second term in Eq. ( |l0l| ) is 
much more involved as a further four-dimensional inte- 
gration and up to hundreds of terms are implied in the 
sum for each hadron species. This troublesome problem 
has been dealt with in refs. jll],[l^] and worked out by 
an analytical calculation of the kernels Kxjk for two and 
three body decays and performing integrations either ana- 
lytically (whenever possible) or numerically for a reduced 



set of hadrons feeding the examined j th hadron species. So 
far, this method has not allowed in practice an exhaustive 
computation of the resonance decay effect onto observed 
spectra. Therefore, we have devised a new method based 
on a mixing of Monte-Carlo and analytical calculations 
that is described in the following. 

Let us first c onsid er one term in the sum in the right 
hand side of Eq. (101); the Boltzmann limit is appropriate 
for all k's as pions, the lightest hadrons, are certainly not 
involved and one can write: 



/ dn,- \ fe ^ 



\dp T 
x Z(T 1 Q_^V) 
Z {T, Q, V) 



d t K Tjk (pT,t) TTT^ c 



(2tt) 3 



0a T -Pj 



(107) 



The above expression is in fact the transverse momentum 
spectrum of hadron j stemming from the decays of hadron 
k emitted from a cluster at temperature T boosted along 
x axis. Instead of trying to calculate all kernels Kxjk, our 
starting point is to express this spectrum as a function of 
the same spectrum in the EGC's rest frame. Let us then 
consider the invariant spectrum: 



/ dn 



dn 
dV 



k lf{p)) 



(108) 



where the starred quantities pertain to the cluster's rest 
frame, p = (e, p) is the four-momentum and p* is linked 
to p through the Lorentz transformation Li(ut) ■ The 
transverse momentum spectrum can be obtained by inte- 
grating the cylindrical coordinates p z and (p, namely: 



/ dn \ 
\ dpx ' 



dPz 



(109) 



The rest-frame spectrum ^dn/d 3 p*^ , expressed as 

a function of p*, has the remarkable feature of being 
isotropic. In fact, no dependence on direction is possible 
for the momentum of a hadron (J) produced by the decay 
of a hadron (fc) whose primary spectrum is supposed to 
be, according to the prediction of the statistical model, 
isotropic. Therefore: 



d 3 p* 



k^j 
(P* 



afe) (p* 

47rp* 2 



&) (O 

47rp*e* 



(110) 



By using Eq. (11C), Eq. (10S) can be written as an integral 
transform: 



/ dn \ _ 
\ dpT I 

= I dx f + °d Pz ["dep ^ ( ^) k lk 6(x - e*( P )) 

J — OO 



_ 7r 47rp*e \ de 



(in) 
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The energy in the cluster rest frame can be calculated by 
using the formula for a Lorentz transformation along x 
axis: 



e*(p) = 7t£ — PtItPt cos ip = y 1 + Uj,e — utPt cos ip 

The integration on the variables p z and <p can be per- 
formed analytically and the result is (see Appendix C): 



dp T , 



Ap 



.2 _ JO Xd P* 



1 + TTlT 



F o. r ( 113 ) 



27rp* ^/ (z + - z mill )(z mal + 1) V 2 ' 
where F is the elliptic integral of the first kind and: 



c)(Zn 



^max 
z± = 



2m in) (^i: 

= max(l, Z-) 

e* ± utPt 
== ( 

1 + mj- 



min(l, z_) 



(114) 



What have we g aine d by expressing the secondary spec- 
trum as in Eq. ( |113| )? In fact, we have unfolded it s dc - 
pendence on ut in an explicit form, unlike in Eq. ( |107| ), 
and this has a strong impact on the problem overall, since 
we have now to calculate only the momentum spectrum 

dn/dp* ) in the cluster rest frame, whose shape de- 
pends only on T. Instead of doing this analytically, which 
is much too complicated for three or more body decays, 
we have adopted a Monte-Carlo procedure which may be 
summarized as follows: 



states quoted in JETSET program tables jl9[ pre- 
dicted by HQET 

2. the obtained momentum spectra of the various hadrons 
emerging from the decay chain initiated by the hadron 
k have been stored in 200 bins histograms ranging from 
to 3 GeV for hadronic collisions and in 333 bins his- 
tograms ranging from to 5 GeV for c + c~collisions, 
owing to the hardness of heavy flavoured hadron de- 
cay products. Thereafter, they have been smoothed 
according to a quadratic interpolation procedure [ pl| , 
normalized to 1 and stored into readable files 

3. the procedure described in 1) and 2) has been repeated 
for temperatures ranging from 140 to 190 MeV in steps 
of 1 MeV. 

Overall, for hadronic collisions, the Monte-Carlo calcu- 
lation has involved the simulation of 200,000 decays of 144 
light flavoured hadronic states for 51 temperature steps, 
yielding a grand-total of 59,721 spectra and w 1.2 10 7 
stored real numbers. For e + e~ collisions, also the decays 
of 83 heavy flavoured hadron decays have been simulated 
yielding a grand-total of 336,447 spectra and « 1.1 10 8 
stored real numbers. 

If i>kj is the momentum spectrum for the decay 
j, normalized to 1, and fkj is the mean number of hadrons 
j produced by the decay chain started by the primary 
hadr on fc, the final f ormu la obtained by adding primary 
(106) and secondary ( 113 ) contributions reads: 



\ dp- 



hi t Pt 



Ki 



1 + uip tut 



vT (utPt\ \ - 4(n fc )P rim / fcj p T 
Xl °l^7-J+^ 

k wl + My rriT 
x ^(p*)F(7r/2,r) 



1. at a fixed temperature, for each light flavoured hadron 
k with mass < 1.8 GeV, we have simulated 200,000 
decays distributed among the known decays modes ac- 
cording to the known branching ratios. If the decay 
products of k were unstable particles, the decay chain 
has been continued until 'stable' states (where 'stable' 
is an experiment-dependent definition) were reached. 
The kinematic distributions of hadronic decays have 
been calculated according to relativistic phase space 
p9[. Resonances with a width r > 1 MeV have been 
given a distributed mass according to a relativistic 
Breit-Wigner (see discussion about Eq. (106)). The 
considered hadron species, their masses and widths, 
their decay modes and relevant branching ratios have 
been taken from ref. pQ|j , For the special case of 

e + e~collisions, where heavy quarks are abundantly 5 Data analysis 



o ' 27rp* yj (z + - z min )(z max + 1) 



(115) 



for all hadrons except pions, for whi ch B ose-Einstein 
statistics cannot be neglected and Eq. fll05p for the pri- 
mary spectrum must be used. The (n) prim 's are the pri- 
mary multiplicities. The fkj have been computed on the 
basis of experimentally k now n branching ratios. The mo- 
mentum integral in Eq. (115) has been computed numer- 
ically by taking into account that Vjk are stepwise func- 
tions vanishing for p* > 3 GeV (5 GeV for e + e~ collisions). 
This truncation introduces a further numerical error in the 
computation which is found to be negligible. 



produced, also heavy flavoured hadron decays have 
been simulated by using known branching ratios [20| 
along with the predictions of the Lund model [19| 
for unmeasured channels. The set of heavy flavoured 
states included all measured ones fl2C(l and all other 



The main goal of data analysis is the assessment of the 
consistency of experimental spectra with the prediction of 
the statistical hadronisation model. This has a twofold im- 
plication: firstly, spectrum shapes of different hadrons at a 
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Table 1. Parameters of the statistical hadronisation model fitted with average multiplicities in hadronic and ee~ collisions. 
Also quoted the Wroblewski factor Ag = 2(ss)/((uu) + (dd)) calculated with the newly produced valence quark pairs from fitted 
primary hadron multiplicities. The very large value of \ 2 in PP collisions is owing to the lack of systematic errors for the 
measurements quoted in the main reference |23J] (see also Table 2). 



Collision 


y/s (GeV) 


T (MeV) 


VT 3 


<ss> 


X 2 /dof 


As 


K+p 


11.5 


176.9 ± 2.6 


5.85 ± 0.39 


0.347 ± 0.020 


68.0/14 


0.203±0.012 


K+p 


21.7 


175.8 ± 5.6 


8.5 ± 1.1 


0.578 ± 0.056 


38.0/9 


0.227±0.028 


7T + p 


21.7 


170.5 ± 5.2 


10.8 ± 1.2 


0.734 ± 0.049 


39.7/7 


0.282±0.032 


PP 


27.4 


162.4 ± 1.6 


14.24 ± 0.66 


0.653 ± 0.017 


315.2/29 


0.200±0.005 


Collision 


(GeV) 


T (MeV) 


VT 3 


7s 


X 2 /doi 


As 


e+e~ 


14 


167.4 ± 6.5 


9.7 ± 1.6 


0.795 ± 0.088 


1.5/3 


0.243±0.036 


e+e~ 


22 


172.5 ± 6.7 


10.6 ± 2.2 


0.767 ± 0.094 


1.0/3 


0.263±0.042 


e+e~ 


29 


159.0 ± 2.6 


17.3 ± 1.5 


0.710 ± 0.047 


29.3/12 


0.228±0.015 


e+e~ 


35 


158.7 ± 3.4 


17.5 ± 1.8 


0.746 ± 0.040 


8.8/7 


0.242±0.017 


e+e~ 


43 


162.5 ± 8.1 


16.2 ± 3.1 


0.768 ± 0.065 


3.0/3 


0.261±0.030 



given centre-of-mass energy should be described by essen- 
tially the same parameters T and ut, as long as v3p -C 1, 
as discussed before; secondly, the best-fit T's should be in 
agreement with the temperature fitted with hadron multi- 
plicities. The latter statement holds provided that hadrons 
do not undergo significant elastic rescattering after their 
formation (or chemical freeze-out), a generally accepted 
belief in elementary collision. In view of this objective, 
a natural requirement for the data set for a given colli- 
sion and centre-of-mass energy, is the existence of a con- 
siderably large sample of both measured transverse mo- 
mentum spectra and integrated multiplicities of different 
hadron species. Furthermore, centre-of-mass energy must 
be high enough to allow the use of a canonical formalism, 
what is expected to occur above roughly y/s w 10 GeV 
(see discussion in Sect. 1). Essentially, four collision sys- 
tems fulfilling these requirementshave been found: K+p at 
= 11.5 and V« = 21.7 GeV, 7r+p at y/1 = 21.7 GeV 
and pp at y/s — 21 A GeV. Though the covered energy 
range is not large, these points should allow to settle both 
of the previous issues. 



strange quarks with respect to u, d quarks |]38f , the extra 
strangeness suppression has been parametrised differently 
here. Indeed, the number of s+s quarks has been con- 
sidered as an additional charge Ns to be conserved into 
final hadrons along with electric charge Q, baryon number 
iV and strangeness S. This means that quantum number 
vectors actually have four components for light flavoured 
hadrons: 



Q = (Q,N,S,N S ) 

q j = (Q j ,N j) S j ,N Sj ) (116) 

Unlike Q, N and S, the number of strange quarks to be 
shared among the primary hadrons is not fixed by the 
initial conditions and may fluctuate. Hence, it is assumed 
that the ss pairs are independently produced from the vac- 
uum with fluctuations governed by Poisson distribution. 
The mean number of ss pairs, (ss), has been taken as a 
free parameter replacing 75. Therefore, the actual formula 
for average primary multiplicities reads: 



5.1 Fit to average multiplicities 

The first step of the analysis is a fit to measured aver- 
age multiplicities of the various hadron species. The fit 
procedure is very much alike that in refs. though 
with some significant improvement. The used formula for 
the primary multiplicities is the integral of the spectrum 
( |105| ), in which only the first term of the series is re- 
tained for all hadrons except pions. Hadron multiplici- 
ties can be reproduced only if a further suppression of 
hadrons with valence strange quarks is introduced. This 
extra strangeness suppression has been implemented in 
refs. |])Q)|5| by means of a phenomenological parameter 
75 < 1 reducing the average primary multiplicities by 
7g s , ns being the number of valence strange quarks, with 
respect to the full equilibrium values. However, in view 
of the constancy of the ratio between newly produced 



(n j ) primary 



VT(2J j + 1) 
2V 2 



M e-^(ss) 



£ 

K=0 



K 



K\ 



n=l 



7) 



with Q and like in Eq. (116) and N s = 2K + N S where 



Ng is the number of valence strange quarks in the col- 
liding particles. The numerical computation of partition 
func tions with four fixed quantum numbers, such as in 
Eq. flll7| ), requires the implementation of a quite involved 
numerical algorithm (described in Appendix D) to keep 
it within practical times. Notwithstanding, it becomes ex- 
tremely slow even for relatively small number of ss p airs, 
say 5. Therefore, the sum over ss pairs in Eq. (117) has 
been truncated to min[3, (ss) + 3\/(ss)] and the Poisson 
distribution renormalized accordingly. As (ss) turns out 
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Table 2. Particle multiplicities in hadronic collisions. 



Particle 



Me 



Fitted 



Reference 



K+p y/S : 

—^n 

7T + 

K+ 
K~ 
P° 
K*° 
K*° 
K* + 
K*~ 
f 2 (1270) 
K*° 
.+ 



11.5 GeV 



2.50 ± 0.12 
2.860 ± 0.087 
1.960 ± 0.025 
0.685 ± 0.062 
0.087 ± 0.025 
0.254 ± 0.020 
0.249 ± 0.023 
0.0685 ± 0.019 
0.253 ± 0.019 
0.0180 ± 0.0093 
0.0734 ± 0.0140 
0.0616 ± 0.0068 
0.0653 ± 0.016 
0.0678 ± 0.0031 
0.0404 ± 0.0025 
0.0124 ± 0.0037 
0.00436 ± 0.00310 



2.724 

2.802 

1.973 

0.6563 

0.09420 

0.3468 

0.2150 

0.03225 

0.2749 

0.02757 

0.05044 

0.02266 

0.03292 

0.06828 

0.04266 

0.01089 

0.005099 



9* ] 

mr b 
S j ab 

£]<• 

S] b 

K] ab 

S*l b 

5 ] ab 

B] ab 

2] ab 

S3] bc 

£] b <= 



K+p \/i ■■ 



charged 
T 



K* + 
f 2 (1270) 



+ + 



-- 21.7 GeV 

8.21 ± 0.12 
2.93 ± 0.47 
0.390 ± 0.045 
0.308 ± 0.037 
0.290 ± 0.028 
0.142 ± 0.027 
0.319 ± 0.042 
0.095 ± 0.017 
0.092 ± 0.008 
0.061 ± 0.009 
0.214 ± 0.034 
0.021 ± 0.009 



8.228 

3.717 

0.4166 

0.4820 

0.2486 

0.06152 

0.2980 

0.06931 

0.1009 

0.06688 

0.1596 

0.01519 



3! ] 



rr+p ^/S = 21.7 GeV 



charged 
77- 

? 

K .o 
R .o 



A + + 



8.31 ± 0.10 
3.092 ± 0.47 
0.236 ± 0.025 
0.460 ± 0.020 
0.129 ± 0.014 
0.148 ± 0.014 
0.103 ± 0.008 
0.032 ± 0.005 
0.221 ± 0.025 
0.015 ± 0.007 



8.434 

4.003 

0.2673 

0.5322 

0.09909 

0.08879 

0.1195 

0.03259 

0.1707 

0.01984 



PP : 
—Ji 



27.4 GeV 



K+ 

K — 
ft' 

K s 

V 

p° 
p+ 

p 

K* + 
K*~ 

f (980) 
4> 

f 2 (1270) 
p 3 (1690) 



/1(1520) 



3.87 ± 0.12 
4.10 ± 0.11 
3.34 ± 0.08 
0.331 ± 0.016 
0.224 ± 0.011 
0.225 ± 0.014 
0.30 ± 0.02 
0.384 ± 0.018 
0.552 ± 0.082 
0.354 ± 0.058 
0.390 ± 0.024 
0.120 ± 0.021 
0.0902 ± 0.016 
0.132 ± 0.016 
0.0875 ± 0.012 
0.0226 ± 0.0079 
0.0189 ± 0.0018 
0.0921 ± 0.012 
0.078 ± 0.049 
1.200 ± 0.097 
0.063 ± 0.0020 
0.1230 ± 0.0062 
0.0155 ± 0.0034 
0.0479 ± 0.015 
0.0128 ± 0.0061 
0.218 ± 0.003 
0.1410 ± 0.0079 
0.0128 ± 0.0049 
0.0335 ± 0.0098 
0.0204 ± 0.0024 
0.0101 ± 0.0018 
0.0171 ± 0.003 



4.594 

4.479 

3.612 

0.3085 

0.1852 

0.2377 

0.4046 

0.5830 

0.6236 

0.4698 

0.4798 

0.09458 

0.06278 

0.1080 

0.05710 

0.03876 

0.02401 

0.06623 

0.009045 

1.054 

0.05277 

0.1461 

0.01669 

0.04369 

0.03252 

0.2514 

0.2057 

0.009645 

0.01426 

0.02060 

0.01396 

0.01054 




a - Only statistical error 
b - The multiplicity has beei 
by the inclastip - Iross sectior. 
collaboration [bd] 
c - Systematic error maybe i 
d - The multiplicity has bee: 
by the inelastic cross sectior 
and cr in = 32.80 mb for pp c 
multiplicities and cross sccti 



culated by dividi 
— 16.07 ± 0.11 



■cd by the sail 



ncludcd in the quoted experimental error 
leulatcd by dividing the quoted cross sccti* 
L = 17.53 mb for K+p, <r in = 20.71 mb for 
sions, inferred from combined quotations ol 
of other particles. 



to be < 1 in all examined collisions (see Table 1), the 
truncation is essentially harmless. 

Once primary average multiplicities of all (light- 
flavoured) hadron species up to a mass of 1.8 GeV are 
calculated for a given set of free parameters T, V and 
(ss), final multiplicities to be compared with experimen- 
tal data are calculated by using the known decay modes 
and branching ratios. The hadronic decay chain is con- 
tinued until particles considered stable by the experi- 
ments are reached. For the presently examined hadronic 
collisions, measured at fixed target, all weakly decaying 
light flavoured hadrons are considered as stable. The var- 
ious hadron species, their masses and widths, their decay 
modes and relevant branching ratios have been taken from 
ref. f§. 

The free parameters T, V and (ss) are determined by 
means of a x 2 fit to measured average multiplicities. For 
each experiment, the most recent measurements have been 
considered. Multiple measurements from different exper- 
iments have been averaged according to the PDG's [^0] 
method with error rescaling in case of discrepancy. The 
minimisation of the \ 2 is performed in two steps in order 
to also take into account the uncertainties on input param- 
eters such as hadron masses, widths and branching ratios, 
according to the procedure described in ref. Q , which is 
summarized hereafter. Firstly a \ 2 with only experimen- 
tal errors is minimised and preliminary best-fit model pa- 
rameters are determined. Then, keeping the preliminarly 
fitted parameters fixed, the variations Zin^ thco of the mul- 
tiplicities corresponding to the variations of the I th input 
parameter by one standard deviation are calculated. Such 
variations are considered as additional systematic uncer- 
tainties on the multiplicities and the following covariance 
matrix is formed: 



r-sys 



= 5> 



(118) 



to be added to the experimental covariance matrix C cxp . 
Finally a new y 2 is minimised with covariance matrix 
C exp + C sys from which the best-fit estimates of the param- 
eters and their errors are obtained. Actually more than 
130 among the most relevant or poorly known input pa- 
rameters have been varied. This fit technique upgrades 
that used in a previous analysis of multiplicities in c + c~, 
pp and pp collisions ||,(|,[|] in that also the off-diagonal 
elements of C sys are included. 

The results of the fit are shown in Tables 1, 2. The 
\ 2, s are not as small as it would have been expected in 
a statistically consistent analysis. Nevertheless, it should 
be taken into account that many measurements only in- 
clude statistical errors and that the large \ 2 value is often 
owing to a single large deviation. Moreover, it should be 
emphasized once more that the assumptions the existence 
of an EGC relies on, may not be fully realised. 

The found values of temperatures and Wroblewski pa- 
rameters As are very close to those found in previous 
analyses j|,[|,|38|; the difference between old and 
present T value in pp collisions at yfs = 27 'A GeV (and 
e + e~collisions as well, see Subsect. 5.3) is mainly owing 
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Fig. 2. x 2 contour lines in the T — ut plane for the fits to iden- 
tified particle transverse momentum spectra in pp collisions at 
y/s — 27 A GeV. The contours range from Xmin + 1 (1 (?) to 
Xmin + 81 (9a's). 



Fig. 3. x' 2 contour lines in the T — ut plane for the fits to 
identified particle transverse momentum spectra in K + p colli- 
sions at yfs — 21.7 GeV. The contours range from Xmin + 1 (1 
<r) to Xmin + 25 (5 ff's). 



to the new parametrisation ((ss) instead of 75), to the up- 
graded fitting procedure taking into account correlations 
in the systematic errors, to the extension of hadron mass 
spectrum cut-off (from 1.7 to 1.8 GeV) and the use of an 
updated set of hadronic data [|(| . A noteworthy feature is 
the sizeable increase of temperature as centre-of-mass en- 
ergy decreases (see Table 1). Due to uncertainties relevant 
to the assumed physical picture and the lack of systematic 
errors in several measurements, it cannot be established, 
for the present, whether this increase is a genuine physical 
effect rather than a numerical artefact in the fit. However, 
it should be remarked that such an effect has been been 
advocated as a possible explanation of the constancy of 
p/7r ratio 



5.2 Fit to transverse momentum spectra 

The mea sure d transverse momentum spectra are fitted to 
the Eq. (115) for all hadron species except pions; for the 
latter, a special treatment was necessary, which is dis- 
cussed later. A peculiar feature of the statistical model is 
the strong relation between multiplicities and spectra as 
both have a chief depe ndence on the same parameter T. 
However, as Eq. (115) shows, final hadron spectra have 
both a direct and an indirect dependence on tempera- 
ture through primary multiplicities. Since we want to test 
the consistency between the chemical temperature T c h ob- 
tained from fits to multiplicities and that - denoted as T - 
obtained from spectrum shape analysis, we have fixe d the 



primary multiplicities (nfc) prim and (nj) pnm in Eq. (115) 
to the values obtained from the independent multiplicity 



fits as shown in Tables 1 and 2. This position also al- 
lows to disentangle the dependence on T of the shapes of 
all spectrum components from that of their overall nor- 
malization. As far as pions are concerned, Bose- Einstein 
statistics cannot be neg lecte d and the first term on the 
right hand side of Eq. (115) must b e re placed with the 
series on the right hand side of Eq. (105), actually trun- 
cated at n — 5. This is quite a special case compared with 
all other hadrons because chemical factors also affect the 
shape of the spectrum besides its normalization. There- 
fore, both integrated primary multiplicity and chemical 
factors of pions have been fixed to those calculated with 
the parameters V,T c h, (ss) of the multiplicity fit. 

Each measured spectrum is fitted with three free pa- 
rameters: T, ut and an overall normalization factor A 



multiplying the right hand side of Eq. (115). In principle, 
this factor would not be needed if the average multiplici- 
ties (rife) and (rij) were in perfect agreement with experi- 
mental measurements. However, this is never the case as 
more or less pronounced deviations from measured values 
usually occur and, also, transverse momentum spectrum 
and multiplicity may have been measured in different ex- 
periments. Therefore, the normalization parameter A has 
to be introduced in order to ultimately decouple the de- 
pendence of spectrum shape on T and ut from that of its 
integral. 

The fit consists of a minimisation of a x 2 for each mea- 
sured hadron j: 



M 



{Ariji thoo - Ai 



ji cxpj 



»=1 



®i exp + ® i sys 



(119) 
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Fig. 4. \ 2 contour lines in the T — ut plane for the fits to 
identified particle transverse momentum spectra in 7r + p colli- 
sions at y/s = 21.7 GeV. The contours range from Xmin + 1(1 
a) to Xmin + 25 (5 cr's). 



Fig. 5. Confidence level contour lines relevant to the x 2 test 
for spectra in pp collisions at y/s = 27.4 GeV. The contours 
range from a confidence level of 31.7% (1 a) to 5.7 10 (5 cr's) 
in steps of 1 a, except K~ for which the innermost contour- 
corresponds to a confidence level of 21.1%. 



where the sum runs over the M bins of the experimental 
spectrum and Ariji are the integrals over the i th bin of 
the measured spectrum: 



Ariji = 



'bin 



(120) 



The exponent k stands for the different variables used for 
the spectra, essentially pr or p\. As far as the experimen- 
tal value is concerned, the above integral is simply the 
product of the quoted differential spectrum times the bin 
width. Besides the quoted ex perim ental error a 



have also included in the x ( 119 ) the systematic error 



l cxpi 



we 



(added in quadrature) due to the uncertainty on masses, 
widths and branching ratios of involved hadrons. This is 
calculated by using the estimated systematic uncertainties 
Sj on the total multiplicity with the method described in 
the previous subsection: 



correlations between uncertainties on p rima ry multiplici- 
ties (pointed out in Subsect. 5.1) in Eq. (115) are neglected 
and this should lead to somewhat optimistic x 2 values. 
This effect can be seen in the residuals distributions of 
fits to pion spectra in Fig. ^ which show some coherent 
structure indicating a bin-to-bin correlation. 

The x 2 minimisation yields preliminary T, ut and A 
values for each hadron. Since both T and ut parameters 
determine the slope of the spectrum, they are expected to 
be correlated variables in the fit. Hence, in order to check 
the fit outcome quality and assess the degree of correlation 
between T and ut, X 2 contour plots in the T—ut plane are 
calculated by keeping A fixed equal to its preliminary best- 
fit value. It must be remarked that best-fit A is almost 
entirely indepe nden t of T and ut as the full spectrum 
integral in Eq. (115) depend only on the multiplicities (rij) 
and (nk), which are fixed in the spectrum fit: 



(121) 



This error is assumed to contribute to the errors of all bins 
independently and proportionally to the spectrum height: 



(Arij 



'% sys 



(122) 



so that the relative systematic error on the spectrum 
height in each bin is constant. With such an assumption, 
the correlations between different bins stemming from the 



(nj) - A [fa)*** + (n k )f kj ] 



(123) 



Indeed, a residual correlation between A and T — ut is 
brought about by the unmeasured portion of the spectrum 
at high pt which, however, gives a very small contribution 
to the overall integral. Thus, the quasi-independence of A 
on T and Ut makes the pattern of local minima of the x 2 
in the T — ut plane an (almost) absolute one. 

Most x 2 contour plots exhibit a pattern of several shal- 
low local minima along a valley in the T — ut plane (see 
Figs. ||, ||, |) demostrating a strong anticorrelation be- 
tween T and ut in the fit. This feature makes the defi- 
nition of a best fit quite difficult; there are indeed many 



22 



F. Becattini, G. Passaleva: Statistical hadronisation model. 



' 0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 

; 0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 

0. 




Log p = -5.9 



i i i i I i i i i I i i i i I i i i i I : i i i i I i i i i I i i i i I i i i i I i i i i I i i ' 

14 0.145 0.15 0.155 0.16 0.165 0.17 0.175 0.18 0.185 



T (GeV) 



' 0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 

; 0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



i- ' : ' • .•/••>••. 

1 K + p \ ; s = 21.7 GeV 

L I, X, 2 = 107/90 

: 1 1 1 1 1 






L K + p Vs = 21.7 GeV N^.; 
~ ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ 


Logp = -l. 

i—i—i— Li 


8 



Fig. 6. Global x 2 contour lines (t op) a nd confidence level con- 
tours for the test variable p in Eq. (124) (bottom) in the T — ut 
plane for the fits to transverse momentum spectra in pp colli- 
sions at yfs = 27 A GeV. The x 2 contours range from Xmin + 1 lisions at yfs = 21.7 GeV. The % 2 contours range from Xmin + 1 



0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.175 0.18 0.185 

T (GeV) 

Fig. 7. Global x 2 contour lines (to p) a nd confidence level con- 
tours for the test variable p in Eq. (124) (bottom) in the T — ut 
plane for the fits to transverse momentum spectra in K + p col- 



(1 a) to Xmin + 81 (9 <r's) whereas the confidence level contours 
range from 10 -6 to 10~ 14 . Also shown the absolute minimum 
(top) and maximum (bottom) and the vertical 1 a tempera- 
ture band determined by the fit to multiplicities. The local x 2 
minimum chosen as starting point of single particle fits is the 
small one beside the vertical band on the left hand side. 



(1 a) to Xmin + 25 (5 it's) whereas the confidence level contours 
range from 10 -2 to 10 -10 . Also shown the absolute minimum 
(top) and maximum (bottom) and the vertical 1 a tempera- 
ture band determined by the fit to multiplicities. The local x 2 
minimum chosen as starting point of single particle fits is the 
absolute minimum. 



local minima with sufficiently low \ 2 or ; equivalently, pro- 
viding a satisfactory confidence level for the x 2 test (see 
Fig. |). 

In order to better define the best-fit values, the search 
of local minima is performed in the T — ut plane for 
the sum of all x|' s at a given centre-of-mass energy, 
thereby assuming the universality of the slope parame- 
ters T and ut for different hadron species. Due to the 
quasi-independence between normalization parameter A 
and (T,ux), this procedure amounts to fit the spectra of 
all measured hadrons at the same time, by using one T, 
one ut and different normalization parameters A\, A2, ■ ■ ■ 
for each hadron species. A suitable minimum for ■ x 2 
can be determined by intersecting the valley of local min- 
ima with the 1(7 band of the temperature T c h extracted 
from the fit to average multiplicities for the same colliding 
system. This is shown in Figs. [| fj], || for pp collisions at 
^fs = 27 .4 GeV, K+p and 7r+p collisions at </s = 21.7 
GeV. It should also be emphasized that for K + p and 7r + p 
collision the local minimum located within the T c ^ band 
coincides with the absolute minimum and, moreover, the 
local minima have apparently moved towards higher T 
with respect to pp collisions, in agreement with the same 
observed trend of T c h (see Table 1); this is a desired indi- 
cation of a correlation between T and T c h- In Figs. ^, 0, 



|3| are also shown the confidence level contour lines for the 
test variable p Q : 

rx 2 3 (T,u T ) 

P(x, 2 (T,UT),n J )= / fix^n^dx 2 
Jo 

where f(x 2 , n j) is the x 2 —distribution for rij degrees of 
freedom. This variable can be used as a suitable test for 
the hypothesis of consistency, with one T and one ut 
value, among the spectrum slopes of different species, in 
that it combines symmetrically the confidence levels of all 
particles. Conversely, the sum of x 2 ' s > albeit the most ap- 
propriate quantity to estimate the assumed common T and 
Ut, cannot be used as a consistency test variable because 
it favours particle spectra with more, and more accurately, 
measured points. The maximum confidence levels for the 
test variable p (whose location does not need to be the 
same as for the global x 2 ) in Figs. |[ @, || are quite poor. 
However, it should be reminded that the used fit model, 
for hadron multiplicities as well as for transverse momen- 
tum spectra, rests on several additional assumptions and 
approximations, besides the basic postulate of local statis- 
tical equilibrium, essentially those listed at the beginning 
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Fig. 8. Global \ 2 contour lines (t op) a nd confidence level con- 
tours for the test variable p in Eq. (124) (bottom) in the T — ut 
plane for the fits to transverse momentum spectra in 7r + p col- 



lisions at */s = 21.7 GeV. The x contours range from Xmin + 1 
(1 a) to Xmin + 25 (5 cr's) whereas the confidence level contours 
range from 10 _J to 10~ 13 . Also shown the absolute minimum 
(top) and maximum (bottom) and the vertical 1 a tempera- 
ture band determined by the fit to multiplicities. The local x 2 
minimum chosen as starting point of single particle fits is the 
absolute minimum. 



Table 3. Best-fit temperatures and mean transverse four- 
velocities ut for each identified hadron transverse momentum 
spectrum in hadronic collisions. 



Particle T (MeV) u T 



X 2 /dof Reference 



K+p S~s = 21.7 GeV 



7T 

k° s 

p° 

K*° 

f 2 (1270) 



181.6 
175.0 
158.0 
174.7 
141.4 
161.7 



0.2101 
0.1485 
0.1631 
0.3155 
0.0218 
0.1553 



19.2/17 

39.3/43 

11.2/10 

3.55/3 

1.23/4 

1.61/3 



7T 1 p 



21.7 GeV 



7T 

P° 

K*° 

K*° 



174.0 
171.4 
162.9 
174.5 
140.0 
140.0 



0.3021 
0.1178 
0.1251 
0.1252 
0.0018 
0.0031 



13.2/17 

25.3/43 

8.00/11 

1.05/2 

4.31/3 

5.36/2 



pp y/s = 27.4 GeV 



7T° 
7T+ 
7T _ 

K+ 
K" 
V 
P° 

f2(1270) 



159.6 
161.0 
165.0 
168.5 
158.3 
174.4 
158.9 
156.0 



0.1948 
0.2140 
0.2179 
0.2324 
0.1960 
0.2260 
0.1971 
0.2004 



22.5/22 

37.5/38 

40.2/38 

26.9/38 
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of Sect. 4, which make theoretical formulae expected to be 
valid only to a certain degree of accuracy. Moreover, the 
use of fixed fitted relative multiplicities of parent hadrons 
in transverse momentum spectrum fits involves a further 
systematic uncertainty, which is very difficult to estimate 
without resorting to a global fit of spectra and multiplic- 
ities at the same time, which is beyond the scope of the 
present work. In view of such considerations, it should not 
be too surprising that rigorous statistical tests do not turn 
out to be fully satisfactory, especially for very accurately 
measured spectra. Nonetheless, it would be desirable to 
have a more direct and intuitive feeling whether the ex- 
amined spectra are in agreement with each other within 
the expected accuracy of the model, whereas the measure 
of consistency by means of a maximum confidence level 
can be too criptic in this regard. This and related issues 
will be discussed in detail in the next section by addressing 
the tot scaling property. 

Once that a local minimum is found through the in- 
tersection of the T c h band with the best-fit valley of the 
global x 2 , ah single spectra are refitted in turn by instruct- 
ing the minimisation algorithm p4| to start from that very 
point with a small initial search step. This allows to deter- 
mine the single-particle x 2 's local minimum closest to the 
selected initial point. The refit results are summarized in 
Table 3. The fitted spectra are shown in Fig. ^ along with 



the calculated cumulative contributions of several classes 
of parent hadrons and residuals distribution. Because of 
the closely spaced local minima in the x 2 graphs, the min- 
imisation algorithm has not been able to produce a reliable 
error estimate. Indeed, in most cases, the found local min- 
ima are so shallow that it is not even possible to define an 
error on the basis of a nearby x 2 = Xmin + 1 closed contour 
line. Thus, we have decided not to quote any error on the 
best-fit parameters in Table 3. 

We have also analysed transverse momentum spectra 
of hadrons in K + p collisions at s/s = 11.5 GeV; the rel- 
evant contour plots of 7r° jlB) and K|g jl6| are shown in 
Fig. [ll]. It can be seen that in both cases the best-fit val- 
ley does not intersect the la T c h band. Therefore, it is 
not possible to find a suitable value of temperature to re- 
produce both integrated yields and transverse momentum 
spectra in this collision system. This failure of can be pos- 
sibly explained by the inadequacy of the canonical frame- 
work at low centre-of-mass energy, where microcanonical 
effects are expected to show up. Particularly, exact total 
transverse momentum conservation should entail a steep- 
ening of the single inclusive transverse momentum spectra 
(a high pt suppression analogous to chemical canonical 
suppression) , an effect which is in agreement with the ob- 
served discrepancy in Fig. [ll| However, this explanation 
is still to be thoroughly tested by carrying out a detailed 
microcanonical calculation. 
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Fig. 9. Transverse momentum spectra of identified particles in pp, K + p and 7r + p collisions at y/s — 27. 4, 2f .7 and 21.7 GeV 
respectively. For each particle, the upper plot is dcr/dp^ (mb/GeV 2 ) versus (GeV 2 ) in logarithmic scale, the lower plot is the 
fit residuals distribution in linear scale in standard deviation units. In the spectrum plots, the upper solid lines are the best-fit 
ones while data are shown as full dots with error bars. Also shown as dashed or solid lines the cumulative contributions due 
to secondary decays for six classes of parent hadrons (from bottom to top: strange baryons, delta-like baryons, N-like baryons, 
strange mesons, non-strange charged mesons and completely neutral mesons). 
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Table 4. Particle multiplicities in e + e collisions. 
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Fig. 10. Fraction of primary pions and kaons in pp collisions 
at \fs = 27.4 GeV estimated by using the fitted parameters in 
Table 3. The relative primary component increases as a func- 
tion of p% but does not reach a value higher than 70-80% even 
at pt as high as 6 GeV 2 . 
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Fig. 11. x contour lines in the T — ut plane for the fits to 
identified particle transverse momentum spectra in K + p colli- 
sions at y/s = 11.5 GeV. The contours range from Xmin + 1(1 
Cr) tO Xmin + 25 (5 a's). 
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5.3 Transverse momentum spectra in e + e collisions 

The very definition of transverse momentum spectrum 
considerably changes from hadronic to e + e~ collisions. In 
the latter, the reference projection axis for px is no longer 
the beam line but a suitably denned event or thrust axis, 
which is the best approximation to the direction of the pri- 
marily created qq pair. However, this axis has to be deter- 
mined experimentally on an event by event basis and this 
is generally done by minimising the sum of tranverse pro- 
jections of particle momenta. Hence, the transverse mo- 
mentum of a particle with respect to the event axis is a 
biased variable in that the same quantity has been used 
to determine the axis itself. Only when number of particle 
per event becomes very large, those spurious correlations 
should be negligible. 

Being aware of this difficulty, we have fitted charged 
particle transverse momentum spectra in e + e _ collisions at 
four different centre-of- mass energies: 14, 22, 29, and 35 
GeV ]57| . Unlike in hadronic collisions, transverse spectra 
have not been measured for several identified particles and 
this compelled us to consider only charged tracks. The fit 
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Fig. 12. x 2 contour lines in the T — ut plane for the 
fits to charged particle transverse momentum spectra in 
e + e~ collisions at y/s — 14,22,29 and 35 GeV. The contours 
range from Xmin + 1 (1 c) to Xmin + 25 (5 a's). Also shown 
the absolute minima and the vertical 1 a temperature band 
determined by the fit to multiplicities. 



procedure and the treatment of the data was essentially 
the same as for hadronic collisions except for the prelim- 
inary multiplicity fit which has been performed with the 
old 7s parametrisation B instead of the aforementioned 
(see Subsect. 5.1) new method. Moreover, the hadronic de- 
cay chain has been extended to include weak decay prod- 
ucts of K°j and hyperons in order to match the exper- 
imental definition of final hadrons in e + e~experiments. 
The inclusion of heavy flavoured events has been imple- 
mented according to the procedure described in ref. j| 
and the branching ratios R q — er(e + c~ — ► qq)/<7ha,d have 
been taken as the lowest order ones R q oc Q 2 . Fit results 
are quoted in Tables 1 and 4. 

The x 2 contour plots in the T — ut plane are shown in 
Fig. IJ] along with the relevant ler T c h bands. Whilst there 
is no intersection between best-fit valley and T h band for 
e + e~ collisions at y/s = 14 GeV, as it was found in K + p 
collisions at y/s = 11.5 GeV, fairly consistent local minima 
can be determined up to y/s — 29 GeV with sufficiently 
low ut values. However, already at y/s — 29 GeV and 
chiefly at y/s — 35 GeV, the fitted Ut is maybe too large 
for the assumed zero-order approximation, described in 
Sect. 3, to hold. This is confirmed by the increasingly high 
X 2 values demonstrating the inability to fit experimen- 



tal transverse momentum spectrum to the shape (115). 
Otherwise stated, as the centre-of-mass energy increases, 
transverse momentum spectra of hadrons become more 
and more influenced by the shape (which can no longer be 
parametrised with an average ut only) of the transverse 
four-velocity distribution of hadronising clusters, which is 



related in turn to the shape of pr distribution of gluons 
radiated in the perturbative QCD shower. Conversely, the 
primordial T-dependent hadronisation spectrum becomes 
less and less important; at high energy, its trace is only a 
small sumperimposed noise smearing out original partonic 
distributions. 



6 Discussion and m T scaling 

As has been discussed at the beginning of previous sec- 
tion, two main predictions of the statistical hadronisation 
model have been provided concerning transverse momen- 
tum spectra: the parameters describing the shape of dif- 
ferent species should be approximately the same and the 
extracted temperature should be in agreement with that 
obtained from fits to average multiplicities. 

As far as the latter statement is concerned, a definite 
answer cannot be given because we have in fact used mul- 
tiplicity fits to pin down the temperature in spectrum fits. 
Yet, an independent indication of consistency is provided 
by the correlation between the temperature location of the 
global x 2 minima in K + p and 7r + p collisions at y/s = 21.7 
GeV and the corresponding la band of T ch (see Figs. [?|, 
||) . A clearcut test would be the analysis of spectra with 
ut — 0, a situation that, in view of the present analy- 
sis, is foreseen to occur in a small energy window below 
\/s k 20 GeV. Focussing solely on the most accurately 
measured spectra, namely charged pions, it can be seen 
that ut appreciably increases from y/s = 21.7 GeV to 27.4 
GeV (from 0.12-0.15 to 0.21) whilst T decreases like T ch 
does, not a surprising fact because the T c h band has been 
actually used to locate a suitable best-fit value for T. On 
the other hand, if pion spectra are singled out in all three 
examined collisions, their mimimum x 2 valleys in Figs. |^, 
||, |] nicely overlap and one might be led to argue that, in 
fact, spectra do no show any significant change of shape 
from 21.7 to 27.4 GeV. In the assumed statistical formal- 
ism, this lack of an observable broadening is due to an 
accidental compensation between a decrease in tempera- 
ture and an increase in transverse four- velocity. This effect 
would be essentially related to a decrease in the chemical 
temperature T c h from 21.7 to 27.4 GeV and an indication 
of its physical genuiness would be the apparent correlation 
between the best-fit global T and T c h band in Figs. 0, @. 
However, for the present, such behaviour of temperatures 
and transverse four- velocities and the possible reason of it 
cannot be firmly established, as discussed in Subsect. 5.1 
with regard to chemical temperature. It might well hap- 
pen that, resorting to a more fundamental microcanonical 
description (with hadronisation occurring, for instance, at 
constant cluster energy density M/V = const), a mono- 
tonic evolution of the parameters governing pt and mass 
slopes as a function of centre-of-mass energy is retrieved 
while the compensation effect between T and ut would 
be only the effect of enforcing a canonical approximation. 
In general, this topic needs a more detailed study. 

Turning back to the first issue, whether statistical pa- 
rameters describing spectra are sufficiently universal for 
different hadron species, we have seen in Sect. 5 that 



F. Becattini, G. Passaleva: Statistical hadronisation model. 



27 



° K 1.5 

J 1.4 

| 1.3 

3 1.2 

J U 

5 1 

e 

3 0.9 
0.8 
0.7 

1 14 
5 13 

3 1.2 

J 1.1 

5 1 

e 

3 0.9 
0.8 
0.7 



L f 2 Final hadrons 


'■. . . i 




T=160 MeV 
u T = 0.2 

J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L 


=■ — f 2 

L ...,p° 
: n 

L ....K + 


Primary hadrons 


: . . . i 




T=160 MeV 
u T = 0.2 

. . i .... i .... i .... i . 



0.5 



1.5 



2.5 3 

m T (GeV) 




0.5 1 1.5 2 2.5 

p T 2 (GeV) m T (GeV) 



Fig. 13. Ratio between mr spectra of quoted particles and 
7T° spectrum taken as reference for primary (bottom) and final 
observable hadrons (top) calculated for T — 160 MeV and 
ut = 0.2 in pp collisions. The spectra have been normalized 
so as to have the same value at the threshold rnr = m. 



Fig. 14. Comparison between experimentally measured p% 
and corresponding tut spectra of pions and kaons in pp colli- 
sions at s/s — 27 .4 GeV. The p% spectra have been normalized 
so as to have the same value in the first bin. Lines connecting 
data points are drawn to guide the eye. 



rigorous statistical tests would rule out this possibility, 
though a fully affirmative answer would have indeed been 
unlikely because of the various approximations inherent to 
the theoretical formulae. However, the implementation of 
a single statistical test does not tell us the main sources 
of the discrepancy between model and data and a fur- 
ther investigation is necessary. As far as pp collisions are 
concerned, it can be realised by looking at Table 3 and 
the location of x 2 minimum valleys in Fig. g, |[ that the 
different fitted temperatures and tranverse four-velocities 
are indeed in good agreement with each other with the 
exception of K + whose spectrum is seemingly harder than 
expected and significantly different from K~'s. The T and 
Ut central values extracted for 77, p° and fi are affected by 
large uncertainties (% 2 minimum valleys are much wider, 
see Fig. ||) so they are essentially in agreement with the 
more accurate ones. The situation is slightly worse in the 
other two examined hadronic collisions, K + p and 7r + p at 
y/s = 21.7 GeV. By inspecting Table 3 and Figs. |, f§ sev- 
eral discrepancies can be noticed: 7r°'s spectrum is harder 
than 7r _ 's, Kg's is softer and K*'s much too soft; these 
features are common to both collisions. 

In order to assess the universality as a function of dif- 
ferent hadron species, a very useful property is the so- 
called mT-scaling, a well known feature of statistical and 
thermal models if St = 0, the shape of the primary 

spectrum of any hadron (except for charged pions at very 
low pt, owing to quantum statistics corrections) depends 
on its mass only thr ough the variable itit — a/Pt + m 2 , 
as shown in Eq. ( 115 ). The mr-scaling is broken by a non- 
vanishing transverse velocity and by the hadronic decay 



chain but both are small effects for the ut values found 
in the present analysis, as it can be seen in Fig. [l3| the 
scaling violation is limited to 25% at the top rriT value 
(~ 2.4 GeV) in the analysed pp collisions. The general 
trend of the scaling violation can be roughly summarized 
with a relative softening of m-r spectra of lighter particles 
with respect to heavier particles. 

The approximate mj- scaling apparently holds in pp 
collisions, for the most accurately measured spectra of pi- 
ons and kaons. That this is a non trivial fact, it is proved 
in Fig. |l4|, where the actually measured and published 
[ p3[ p T spectra are compared with the corresponding spec- 
tra transformed in the variable mj. While p T slopes are 
markedly different between pions and kaons, their mx 
slopes are much closer (heavier particle spectra are af- 
fected by too large errors to allow any conclusion). To a 
closer look, also the aforementioned discrepancy between 
K + and K~ and the hardness of K + spectrum show up; 
this is better seen in Fig. |l5| where predicted and measured 
ratios between mj- spectra of K + and K~ with respect to 
7r + 's spectrum are shown. Although the model succeeds 
in predicting a slight difference between the slopes of K + 
and K~ (in the right direction) owing to different reso- 
nance feeding (the colliding system is not isospin sym- 
metric) , the observed ratios with tt + spectrum are steeper 
in shape, particularly for K + , and the difference between 
K + and K~ is larger. It is very difficult to understand the 
reason of this discrepancy; at present, it can attributed to 
deviations from the model scheme (i.e. from the assump- 
tions listed at the beginning of Sect. 4) or to a systematic 



28 



F. Becattini, G. Passaleva: Statistical hadronisation model. 




0.2 0.4 0.6 0. 



1 1.2 1.4 1.6 1.8 2 

m T (GeV) 



Fig. 15. Ratio between mr spectra of kaons and tv + spec- 
trum taken as reference. The experimental ratios (dots with 
error bars) are calculated by using the measured kaon spectra 
and the fitted ir + spectrum (parameters quoted in Table 3), 
which is a very good approximation to the measured one. The 
theoretical ratios are calculated by taking common values for 
the relevant parameters, namely T = 160 MeV and ut = 0.2, 
which are a fair average of values quoted in Table 3. The spec- 
tra have been normalized so as to have the same value at the 
threshold mr = itik- Lines connecting data points are drawn 
to guide the eye. 



error in the estimation of relative abundances of parent 
hadrons. 

A similar behaviour is observed in K + p and 7r + p colli- 
sions: as shown in Figs. [l6j [l7[ the slopes of mr spectra are 
definitely much closer than those of p T spectra. However, 
in agreement with the previously observed contour plots in 
Figs. [| [|and to the values quoted in Table 3, tt° spectrum 
has a remarkably different shape at low mx with respect 
to 7T~'s and K*'s spectra are definitely steeper than ex- 
pected. The softness of K*'s spectrum could be related to 
the presence of a valence strange quark in the initial state 
for K + p collisions; yet it is there in 7r + p collisions too. 
Unfortunately, the spectra in both K + p and 7r + p colli- 
sions have been measured by the same experiment, NA22, 
and this prevents a cross-check on the data to assess the 
genuiness of such similarities. Furthermore, there is little 
overlap between the particle sample in pp and K + p, 7r + p 
with the exception of ir° and 7r~. In this regard, whilst 
in pp collisions very good agreement is found between the 
measured slopes of the three pion species, a considerable 
discrepancy is observed in in both K + p and 7r + p collisions 
with different initial colliding particles. Indeed, this appar- 
ent violation of a nice scaling observed at a centre-of-mass 
energy only 6 GeV higher is quite odd. 
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Fig. 16. Comparison between experimentally measured p% 
and corresponding mr spectra of identified particles in K + p 
collisions at y^s = 21.7 GeV. The p\ spectra have been nor- 
malized so as to have the same value in the first bin. Lines 
connecting data points are drawn to guide the eye. 
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Fig. 17. Comparison between experimentally measured px 
and corresponding mr spectra of identified particles in 7r + p 
collisions at ^fs = 21.7 GeV. The p% spectra have been nor- 
malized so as to have the same value in the first bin. Lines 
connecting data points are drawn to guide the eye. 
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7 Summary and conclusions 



Appendix 



Transverse momentum spectra of identified hadrons in 
several high energy collisions have been compared to the 
predictions of statistical hadronisation model, which has 
been described in detail starting from a microcanonical 
formulation, with emphasis on the assumptions needed 
to introduce the usual concepts of single-cluster tem- 
perature and volume. The distortion effect due to sec- 
ondary hadronic decays has been accurately and exhaus- 
tively taken into account by a newly proposed numerical 
method. A good agreement is found between model and 
data and, on top of that, the temperatures estimated by 
fits to average particle multiplicities and those extracted 
by fitting transverse momentum spectra are fully compat- 
ible with each other. This is an indication in favour of 
one of the key predictions of the statistical hadronisation 
model, namely the existence of a close relationship be- 
tween the production of particles as a function of their 
mass and, for each particle species, the production as a 
function of momentum (measured in the rest frame of the 
cluster they belong to) at the hadronisation. 

The model calculations have been performed by tak- 
ing several assumptions. Firstly, the canonical framework 
has been used, which is expected to fail at low centre- 
of-mass energies where the effect of exact transverse mo- 
mentum conservation must be significant and, in fact, in 
e + e~ collisions at a centre-of-mass energy of 14 GeV and 
K + p collisions at 11.5 GeV, a clear discrepancy between 
data and calculations shows up. Secondly, the approxima- 
tion of very small average transverse velocity of the clus- 
ters has been introduced, which should break down at high 
centre-of-mass energy because of the increasingly large pr 
parton emission (giving rise to clusters with large tran- 
verse velocity). At very high energy, this effect is predom- 
inant in determining the shape of transverse momentum 
spectra whereas hadronisation plays the role of a small su- 
perimposed smearing noise. Indeed, it is found that even at 
moderately high energy e + e _ collisions, at s/s — 35 GeV, 
the approximate formulae are no longer able to reproduce 
the shape of experimental charged particle spectrum and 
the resulting transverse four- velocity is not consistent with 
the initial requirement. In summary, the presently used 
parametrisation of px spectra is found to work well only 
in a limited centre-of-mass energy range (roughly between 
20 and 30 GeV). 

In order to further investigate statistical features of 
hadronisation, the study of px spectrum at lower centre- 
of-mass energy would be greatly interesting but full com- 
plex microcanonical calculations are required. 
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A Proof of the reduction to the EGC 

We want to prove by direct transformation that : 



n 
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[ n e / a4p * q ] s ^ p - ^ <^q. 

_ (A.I) 

is the density of states of an EGC with four-momentum 
P = quantum numbers Q = 53iQ* an d vol- 

ume (in the reference frame where four-momentum is P) 
V = Vi . Let us consider the integral representation of 
Dirac's and Kronecker's S: 



S\P - Z t Pi) 



1 



(2tt)' 
1 



lim , 
77^0 (2tt) ? 



d 4 x jl p -* p *> 



d n 6 jiQ-ZiQ^ t-Zirii Q* 



(A.2) 



where rj = (r]i, . . . , r/ n ) is a set of real positive numbers 
which are introduced in order to regularize infinite sums 
for quantum numbers running from to +oo, suc h as the 
absolute number of s trange quarks. By using Eq. ( A.2 ) we 
can rewrite Eq. (A.l) as: 
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Let us focus on the single i th factor in the above product. 
If one plugs in it the expression (0) of J7j Q and the 9 
function's integral representation: 



1 



2tti J-oo-ie " 



+oo-ie _iaP" 



(A.4) 



the following expression is obtained: 



+oo — ie 



JaP°-iPi-(x-Xi)-iQi-(<t)-<t>i)-Vi-Qi 



(A.5) 



where t he li mit r\ — > is implied. The integration in Pi 
in Eq. (|A.5| ) formally (that is, in a distributional sense) 
yields (27r) 4 5(x° — x° + a)S 3 (x.i — x). On the other hand, 
the sum over the quantum number vectors Qi yields, for 
its generic fc th component: 

4 the integration over the line (— oo — ie, +oo — ie), with s — *• 
in the complex plane for all 's, as well as for x°, is implied. 
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27r5(0 fc - <f>ik) if Q lk = — oo, . . . , +00 
1 



m™o 1 - exp[i(0 ifc - <j> k ) - r] k ] 



if Q lk = 0, . 



(A.6) 

In fact, the second expression in the above equation is also 
equivalent to a 8 distribution. This can be shown through 
a transformation of the integration in fa over the interval 
(—it, tt) onto the unitary circle in the complex plane by 
setting w = exp[— ifa k ]: 



lim 



1 



lim — 

>)fc^0 1 



dw 



- exp[i(fa k - 
exp^fx^c 



- mi 



>0n)] 



. Q—iipk — rik 



(A.7) 



The function F w (see Eq. (||)) is analytic with respect to 
the complex variable w as only positive integer powers of 
w are involved (the gjy's are positive integer numbers by 
assumption) and the pole e~ 1 ^ k ~ r,k lies within the unitary 
circle. Thus, the result of the integral in Eq. (A.7) is sim- 
ply: 



lim 27rexp[F K; 



-i<Pk-rik 



,fan)} 



= 2<7rexp[F(Xi,fai,...,<f)k,...,fa n )] (A. 8) 

which is the same result that would have been obtained 
by integrating with a 5 distribution. Of course, the same 
procedure may be repeated for all components othe r than 
k. Eventually, the integrations in Xi and <p i in (|A.5|) yield: 



lim 

e^O 27T1 



da 



exp[i^(x° — a, x, </>)] 



(A.9) 



This integral can be calculated by closing the z = — ie line 
with a semicircle lying in the upper half of the complex 
plane and taking the limit for the circle radius going to 
infinity. Since the function F is analytic with respect to a 
(see Eq. @), the result is simply: 

exp[F(a; , x, </>)] = exp[F{x, 0)] (A.10) 

which is the final form of expression (|A.5| ). Now the func- 
tion F depends on the particular i th cluster only through 
the volume Vi with a simple linear relationship, namely 
F{x, (p) = Vif(x, cj)) (s ee Eq. (||)). Therefore, the product 
of all the factors ( |A.5|) yields: 



exp^Fi/^,0)] 



(A.ll) 



and this is the fi nal f orm of the expression between square 
brackets in Eq. (A. 3). Hence, the fl function reads: 



fl = 



1 



(27T) 4 +™ 



d 4 x / d^e'^^e^)^ 1 '^ 



QED. 



(A.12) 



B Decomposition of Lorentz transformations 

A general Lorentz transformation A can be uniquely de- 
composed as the product of a transformation like in 
Eq. ( f42| ) and a general rotation jlTj : 

A = LR L = ftsfa) R 2 (0) Cs(0 (B.l) 
where <p G [0, 2tt), 9 G [0, tt] and £ > (see Eq. @). 
Any timelike vector can be uniquely written as L(t) where 
i = (1, 0, , 0) is the time axis unit vector and L like 
in Eqs. (^2,BT) jL7|. Hence, given two timelike vectors 
P' = L'(t) and P = L(t) t here is only one Lorentz trans- 
formation of the type ( 42 |B . 1 ) transforming P' into P, 
namely L" such that L"{t) = LL'-^f). 

In order to prove the alternative decomposition of L in 
Eq. ([79|), it is sufficient to show that the equation: 



ftafoO R 2 (0) C 3 (0(t) = L 3 (r?) R 3 ($ ^(0(4) (B.2) 

with 6 [0, 27r), 77 6 (—00, +00) and £ > 0, establishes a 
biunivocal co rresp ondance between the two sets of coor- 
dinates. Eq. (B^) leads to the system of equations: 



cosh £ cosh 77 = cosh £ 
sinh £ cos <p = sinh £ sin 9 cos 
sinh £ sin = sinh £ sin sin 93 

cosh Q sinh 77 = sinh £ cos 

whose solution is: 

cosh £ = cosh £ cosh 77 



(B.3) 



tan# 



tanh£ 
sinh 7/ 



(B.4) 



and its inverse: 



cosh£ 



(B.5) 
(I? 



cosh 77 = — , 

V 1 + sinh 2 £ sin 2 9 

sinh £ = sinh £ sin 

(f> = (p 

Therefore, if Eq. ( f42| ) is a valid decomposition, Eq. 
holds as well, within the relevant variable domains. 

The associated invariant measure in the variables £, f) 
and <f) can be ob tain ed by calculating the jacobian of the 
transformation (B.4) and taking into account that dL = 
sinh 2 £ sin 9 d£ d9 dip [18| . It is advantegeous to rewrite 
the above measure as dL = sinh £ d cosh £ d cos 9 dip and 
use t he foll owin g transformation formulae derived by Eqs. 

flO and (Pjl) : 



cosh£- = cosh C cosh 77 

cosh C sinh 77 

cos# — 



\J cosh 2 £ cosh 2 77 — 1 
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implying: 



! det J| 



cosh £ sinh ( 



\J cosh 2 C, cosh 2 r\ — 1 



and, consequently: 



dL = — — — - d£d77< 



(B.6) 



(B.7) 



C Calculation of the secondary contribution 
to transverse momentum spectrum 



The integration of Eq. (Ill) is performed by solving the S 
function with respect to the tp azimuthal angle according 
to the well known formula: 



(C.l) 



where </?oi are the zeroes of the / function. In the present 
case (see Eq. (|l 1 2j ) ) : 



f((p) — x — 7t£ — utPt cos (p 



thus: 



cos (fo 



\f'M\ 



7re 



utPt 



(C.2) 



(C.3) 



utpt 



COS z ifio 



There are two angles fulfilling Eq. ( JC.3D in the [0, 2tt) in- 
terval, provided that the absolute value of its right hand 
side is < 1. For these angles, the absolute value of the 
derivatives in the lowe r equation are equal. Hence, the in- 
tegration in Eq. (Ill) in the azimuthal angle yields 2 i f 
Ittc — x\ < utPt and otherwise. Therefore, Eq. (Ill) 
turns into: 



where advantage has been taken of the fact that the inte- 
grand is even in rapidity. The last needed change of vari- 
able is z — coshy, leading to: 



/ dn 
\dpT 



-t-oo 



PT 



dx dz - 

'i 4ttJx 2 



\de* 



[z — z-)6(z + — z) 



(C.6) 



7tTOt \/(z — l)(z + l)(z - z-)(z+ - z) 
where z+ and z_ are defined in Eq. (|114|) and, obviously, 



7t = y 1 + u T . The integral in the z variable is an elliptic 
integral whose solutions depend on the sign of Z— — 1; they 
can be found in ref. l59|: 



+oo 

dz 



(z-z-)6(z+-z) 



y/(z-l)(z+l)(z-Z-)( Z+ -z) 



= Pf-,r 
1) V2' 



(C.7) 



where z m - m , z max and r are quoted in Eq. (114) and F 
is the complete elliptic integral of the first kind. After 

the substitution of x with momentum p* = 



x 2 — m 2 as 



integratio n var iable in Eq. (|C.6| ), the final result is exactly 
as in Eq. (pT3|). 



D Calculation of the partition function with 
fixed number of strange quarks 

In this section we describe the technique used to handle 
the numerical integration of canonical partition functions 
with four fixed quantum numbers, namely electric charge 
Q, baryon number B, strangeness S and num ber of s + s 
quarks N$- This can be written as in Eq. (102) with n = 4 
and A, = 1: 



-oo /*+oo 



PT 



-oo Aire* x 2 — m 



(utpt - hrt - x\) 
u\p 2 T - (7 T e - x) 2 



(C.4) 



The next integration can be done with the change of vari- 
able p z — my sinh y, i.e . by introducing rapidity, so that 
e = tot cosh y and Eq. ( p.4[ ) reads: 



PT 



dn . k^j f+°° f+°° 

- — ) = 4 / da; / dy - 

A PT> Jm 3 Ja ATrJ.r' -nr' 

0(utPt — |7t?71t cosh y ~ x\) 
y u T p T — (7ttot cosh y — x) 2 



\de* 



(C.5) 



Z(Q) = -Lj ["dU cxp [iQ • <\> + F c (4>)} (D.l) 
with Q = (Q, B, S, N s ) and: 



F c (<f>) 



V 



(2tt) 3 ^ 

3 



(2Jj + l)J d 3 p log [1 ± e -ei/T-itU-*\±l 

_ (D-2) 

Straight numerical four-dimensional integration is too 
time consuming to allow reasonably quick multiplicity fits: 
an analytical reduction of the above integral is indeed nec- 
essary. For this purpose, we first integrate out the 4>b vari- 
able associated with baryon number. This can be done 
analitically because of the presence of o nly single-valued 
baryons in the function F c in Eq. ( p.l| ), neglecting the 
corrections due to Fermi statistics which are very small. 
Thus: 
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F c (<f>) = F c {4>Q,(j)s,(j)Ns)r 

where: 



W+e _i * B + W-e 14 



(D.3) 



W ± = Z J ex P [-i(<3j>Q + Sj4> s + Nsj^Ns )] (D.4) 

bar. 

*i = V{2 ^ l) J d 3 P cxp [- £j /T] (D.5) 

It must be stressed that W + ^ Wl because the num- 
ber of strange quarks is the same for both parti cles and 
antiparticles. The integration in <pg in Eq. (D.l) can be 
performed by using the decomposition ( D.3 ) and a series 
is obtained: 



J- / dcj) B e iB ^ B exp [W_|_e _i ^ B +W / _e i ^ B ] 

Z7T 



~ w k+B w k 

k\{k + B)\ 

_ fc "° (D.6) 
Then, the integration in 4>n s in Eq. (D.l) is performed. It 
is advantegeous to set w = exp[— i4>N s ] an d map the inter- 
val [— 7T, 7r) onto the unitary circle in the complex plane: 



27ri 



dw 



w- Ns E(w) exp[F M H] 



(D.7) 



where: 



F M (w) = F c ((/)Q,(j)s, 4>N S ) mesons 



(D.8) 

and Z 1 is the series in Eq. (D.6). Since mesons can contain 
at most two strange quarks, the function Fm{w) can be 
written as the sum of three terms: 

F M (w) = a + (3w + jiv 2 (D.9) 

where a, p and 7 are the sums EjZj exp \—'\Qj4>Q — iSjfis] 
over mesons with 0, 1 and 2 strange quarks respectively. 
As baryons may contain at most th ree s trange quarks, the 
functions W+(w) and W-(w) in (D.4) are third-degree 
polynomials in w with conjugate coefficients Aq,...,A3 
depending on <j>Q and <f>s, i.e.: 



3 

W + {w)=Y,A n w r ' 

n=l 



W-{w) 



(D.10) 



and, consequently, the series E{w) in Eq. (|D.6| ) is an ana- 
lytic function of w. Hence, the integrand in Eq. ( p.7|) has 
one pole of Ng 1 order in w = and the integral (D.7) 
turns out to be: 



(27T) 2 



t>s e 



D N ^E(x) exp [p^Q, <p s )x + l{(j> Q As)x 2 ]\ x = 
N S \ 



(D.12) 



which is to be integrated numerically within its two di- 
mensional domain. 

The next problem is ho w to calculate the derivative in the 
integrand in Eq. ( ]D.12 ). We first notice that the series E 
can be calculated analytic ally for a real argument because 
W+(x) = Wt{x) (see Eq. QD.IOD ): 



W+(x) k+B Wl{x) k 



k=0 



k\(k + B] 



l B (2\W+(x)\)exp [iB arg W+(x)} 

(2W+(x)) B (D.13) 



I B (2\W+(x)\) f „ B 



(2\W + {x)\Y 

where 1 B is the modified Bessel function of order B. The 
derivative of order Ns is then expanded: 



D" s E(x)exp[f3x + -/x 2 ] 



N s 



= 12 ( N h) Dk exp $ x + lx 2 ]B Ns - k E(x) (D.14) 



fc=0 



V k 



It can be shown with a little algebra that: 



D fe exp[/3a; + 7x 2 ] 



-«- M ^ (*-*w-fl' (m5) 



The derivative of E(x) can be further expanded: 



The derivative of the first factor in the above sum can be 
calculated by taking advantage of recurrence relations for 
Bessel function derivatives: 



D 



l B {2\W+{x)\) I B+1 (2\W + (x)\) 



(2\W + (x)\) B (2\W + (x)\) B + 
while for the second factor: 



2B\W + (x)\ 2 (D.17) 



exp [a] 

d Ns 

N s \ dz N * 



exp [j3w + "/w 2 ]E(w) 



w=0 



(D.ll) 



The derivative can be calculated along any direction and 
particularly alon g the real axis. Finally, the partition func- 
tion in Eq. ( D.l ) reads: 



B{2W+{x)) B = B(2W+(x)) B - 1 DW+O) (D.18) 

The two equations above are the starting poin t for the 
numerical calculation of the derivative in Eq. ( D.16| ); a 
recursive algorithm has been implemented that comp utes 
it for x — 0, taking into account that (see Eq. ( D.10|) ): 
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l B (2\W+(0)\) _ l B (2\A Q \) 
(2|W+(0)|)* 

D'|W+(a;)| 2 
B l W+(x) 



(2\Ao\) B 
i 



n=0 



x=0 



11 Ai 
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